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1.    Introduction. 

Driving  point  functions  of  ladder  networks  are  frequently  used  in 
the  analysis  and  synthesis  of  these  networks.   The  input  impedance  and/ 
or  input  admittance  can  be  readily  expressed  in  the  form  of  a  continued 
fraction.   Many  good  texts  have  been  written  which  thoroughly  explore 
the  many  facets  of  the  problem.   Recently  the  application  of  the  theory 
of  continuants  to  ladder  networks  has  been  receiving  renewed  interest, 
as  evidenced  by  a  number  of  papers  in  the  current  literature.   Contin- 
uants, provide  a  convenient,  concise,  and  readily  calculable  form  for 
deriving  and  expressing  results  of  network  analysis. 

It  is  often  of  interest  to  view  the  ladder  network  as  a  two-port 
system.   If  we  impress  a  driving  function  at  the  input  terminals,  we  gen- 
erally want  to  know  what  we  will  get  at  the  output  terminals.   Conversely, 
we  might  wish  a  particular  response  from  a  given  driving  function.   The 
task  is  then  to  find  a  network  which  will  give  that  response.   The  prob- 
lem of  synthesis  is  generally  more  difficult  than  that  of  analysis. 

This  paper  shall  be  concerned  with  the  synthesis  of  RC  ladder  net- 
works of  the  first  Cauer  form  from  a  given  voltage-ratio  transfer  function, 
Continuants  are  used  in  the  analysis  of  general  ladder  networks  to  illus- 
trate their  usefulness  in  expressing  transfer  and  driving  point  functions 
in  concise  and  compact  forms.   A  procedure  is  developed  for  synthesizing 
networks  up  to  nth  order.   It  is  shown  that  there  are  many  sets  of  com- 
ponent values  which  can  be  used  to  build  a  suitable  network.   Two  pro- 
grams have  been  written  and  included  to  illustrate  how  one  can  use  a 
digital  computer  to  great  advantage  in  the  synthesis  problem. 


2.    Continued  Fractions. 

A  continued  fraction  is  an  expression  of  the  form 

+     bi 

a,  + 


1   ,+!i_  (2.1) 

2   a3  + 


which  may  be  written  in  a  more  compact  form  with  dropped  +'s  as 


bl     b2    b3 
a,  + 


1   a2  +  a3  +  a  + (2.2) 

In  this  paper  only  continued  fractions  whose  b  terms  are  unity  will 
be  considered. 

A  terminating  continued  fraction  is  one  which  has  a  finite  number 
of  denominators  as  opposed  to  an  infinite  continued  fraction  which  has 
an  infinite  number  of  denominators.   A  continued  fraction  terminating 
in  a  finite  number  of  denominators  defines  a  convergent  of  the  contin- 
ued fraction.   This  convergent  is  usually  expressed  as  a  ratio  p/q.   In 
the  continued  fraction 

111 
a  +  _i_.  _L^  _L^ (2.3) 

1   a„  +  a3  +  a,  +  . . . . 


Pl   a 

— -  =   1  is  called  the  first  convergent: 

ql 

p2   a    1-    al  a2  +  l 

—  =   1  +  r— ■  « is  called  the  second  convergent: 

q2        a2      a2 


Po  i     i    ai  ao  ao  +  ai  +  a0 

r3  a,     1     1    1   2  3    1    3       , ,  .   . 

—  =   1  +  r  —   =  — ; — is  called  the 

q3  a2  +  a3      a2  a3  +  1 


third  convergent;  and  so  on.   The  nth  convergent  is  given  by 
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£n  _1 _1 1_  . 

q       1     a_  +  a„  +     "a 
n  2     3         n 

The  following  properties  of  continued  fractions  and  convergents, 
taken  from  the  text  by  Wells,   are  submitted  without  proof. 

1)  Any  ordinary  fraction  in  its  lowest  terms  may  be  converted  into  a 

terminating  continued  fraction. 

Pn 

2)  The  difference  between  two  consecutive  convergents  — 

qn 
Pn+1    .       1 


and   is 


Vl        qn  qn+l 

3)  The  even  convergents  are  greater,  and  the  odd  convergents  less  than 
the  fraction  itself. 

4)  Any  convergent  is  nearer  than  the  preceding  convergent  to  the  value 
of  the  fraction  itself. 

3.    Continuants. 

The  theory  of  continued  fractions  is  greatly  simplified  by  means 
of  a  class  of  rational  functions  defined  by  a  sequence  of  linear  equa- 
tions.  These  functions  called  Continuants  were  originally  due  to  Euler. 
A  good  discussion  of  continuants  appears  in  the  text  by  Bartlett.  LIj 

The  rational  integral  function,  p  ,  of  the  n  quantities  a  ,  a_,  a„, 

a   is  defined  by  the  set  of  equations 

n  J 

PO  =  l 

pl  =  al 

p2  =  a2  Pl  +  pQ 

W.  Wells,  Advanced  Course  in  Algebra  (New  York:   D.  C.  Heath 
and  Company,  1904),  pp.  449-456. 
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p3   =  a3  p2  +  Pl  (3a) 


P    =  a   P   i  +  P  n 
rn     n  rn-l    rn-2 

where  ithe  coefficients  a.,  a_,  a„ ,  --a  and  the  first  two  functions, 

1   2   3'    n 

pQ  and  p   are  given.   The  function  pn  thus  defined  is  termed  a 
"Simple  Continuant  of  the  nth  Order",  and  is  denoted  by  k(a1 ,  a~,  a„,=-a  ) 
which  may  be  further  shortened  to  k(a1 ,  a  )  for  the  sake  of  brevity  by 
using  only  the  first  and  last  arguments  of  the  function  in  the  notation. 

Simple  continuants  are  then  defined  by  the  set  of  equations 

pQ  =  k(0)   =   1 


Px  =  k(aL)  =  ax 


(3.2) 


P2  =  k(a!>  a2)   =  a2  Pl  +  P0  =  a2  k^al^  +  k^ 


■  k(ax,  a3)   =  a3  p2  +  px  =  a3  k(a1_,  a2)  +  k(ax> 


Pn  =  k(al>  V   =  an  Pn-1  +  Pn-2  =  \  *<aV    Vl> 
+  k(a1}  an_2) 

Each  succeding  function,  therefore,  has  been  defined  in  terms  of  a 
linear  combination  of  the  two  preceding  functions  with  the  exception  of 
p  and  p1 .   It  is  important  to  note  that  the  value  of  the  zero  order 
continuant,  k(0),  is  unity.   If  the  coefficient  a  has  the  value  of 
zero,  then  the  first  order  continuant  k(a..)  also  has  the  value  of 
zero  and  should  be  written  as  k(a   =  0)  =  0  in  order  to  avoid  confusion. 

12 


In  general,  if  the  first  argument  of  the  function  is  a  ,  the 
continuant  for  n  >  r  shall  be  written  as 

p   =  k(a  ,  a  )   =  a  k(a  ,  a   )  +  k(a  ,  a   ) 
n        r   rt       n    r'   n-1       r'   n-2' 

where  p   is  defined  by  the  equations 

p      ,      =      1      =     k(a    ,    a      n) 
r-1  v  r        r-1 

p        =     a        =     k(a    ,    a   ) 
rr  r  r        r 

=     ar+l   Pr  +  Pr-i      =     k(a-'    *^      =     a~J-i    k(a-'    a-} 

+  k(ar,    a       ) 


(3.3) 


p        =ap1+p0     =     k(a    ,    a   )    =  a     k(a   ,    a      .. ) 
*n  n  *n-l        *n-2  v  r'      n  n  r'      n-17 

+  k(a    ,    a        ) 
r        n-2 

A  second  rational  integral  function,  q  ,  may  be  defined  in  a 
similar  fashion  such  that  q   is  the  same  function  of  a_  as  p   is  of  a-,  ; 


q0  is  the  same  function  of  a_,  a„  as  p»  is  of  a, ,  a  •  q   is  the 
3  2    3     2        1    2   nn 


same 


function  of  a„,  a0 a  as  p   ,  is  of  an  ,  a„ a   ..  :  etc..   These  q's 

2'   3     n    rn-l        1    2    n-1 

are  also  continuants  and  may  be  written  as 
qt      =   k(0)   =   1 


q2      =     k(a2)      =     a2 

q3      =     k(a2,    a3)      =   a3   q2  +  ql      =     a3   k(a2'   +  k^0/ 


(3.4) 


q        =     k(a„,    a)      =     aq1+q„      =     a     k(a„,a      .. ) 
nn  2'      n  n   nn-l        nn-2  n  2      n-1 

+  k(a2,    an_2) 


13 


In  relating  the  continuants  which  have  just  been  defined  to  the 

successive  convergents  of  a  continued  fraction,  the  following  expres 

sions  can  be  written: 
P]_   k(a]_) 

al  =  ^  =  wr 


k(a2)  <3'5> 


11  1  Pn  k(al'    an} 

a,    + 


1        a.  +        a„   +      ...    a  q  k(a        a   ) 

2  3  n  n  2'      n' 

Expansion  of  a  continuant  of  fairly  high  order  by  means  of  the 
fundamental  definitions  may  be  a  long  process.   Several  different 
techniques  for  expanding  continuants  have  been  devised  to  decrease  the 
amount  of  work  required.   Bartlett  presents  a  systematic  technique  based 

on  a  schematic  arrangement  due  to  Hindenburg  (1741-1808).  *■   Parker, 

2 
Peskin,  and  Chirlian  have  shown  a  means  of  representing  a  continuant 

as  a  signal-flow  graph.   Expansion  of  the  continuant  is  accomplished  by 

evaluating  the  signal-flow  graph. 

An  alternate  expression  for  the  simple  continuant  k(a  ,  a  )  is 

1    n 

given  the  Sylvester-Muir  determinant: 


1 

A.  C.  Bartlett,  The  Theory  of  Electrical  Artificial  Lines  and 

Filters  (New  York:   J.  Wiley  and  Sons,  1931),  pp.  45-46. 

2 
S.  R.  Parker,  E.  Peskin,  and  P.  M.  Chirlian,  "Continuants,  Signal 

Flow  Graphs,  and  Ladder  Networks",  Proceedings  of  the  IEEE,  Vol  54, 

March  1966,  pp.  422-423. 
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k(ai'  an}   = 


1 

+1 

0 

0 

0 

0 

.      .      0 

0 

0 

1 

a2 

+1 

0 

0 

0    .      .      , 

.      .      0 

0 

0 

0 

-1 

a3 

+1 

0 

0    .       . 

.      .      0 

0 

0 

0 

0 

-1 

a4 

+1 

0    .      . 

.      .      0 

0 

0 

• 

• 

• 

• 

• 

t        ■             • 

.       .      0 

0 

0 

0 

0 

0 

0 

0 

0    .      . 

.     .  -1 

Vi 

+1 

0 

0 

0 

0 

0 

0    .       . 

.       .      0 

-1 

a 

Expansion  of  the  determinant  yields  the  expansion  of  the  continuant. 

The  total  number  of  terms  in  the  expansion  of  an  nth  order 
continuant  is  given  by  the  expression 


1+(n.U   +  (n-2)  (n-3)  +  (n-3)  (n-4)  (n-5)  + 


(3.6) 


which  can  be  proved  equal  to 

n+1 


n+1 


1 


1  +  Z1"   > 


l-V*  ) 


2n+V~ 


(3.7) 


Results  of  these  expressions  can  be  used  as  checks  on  the  expansions. 
The  continuant  expansions  of 

(^(s)   =  k(Rx,  sCn) 

for  values  of  m  up  to  m=8  are  given  in  Appendix  IV.   These  will  be  used 

in  subsequent  sections  of  this  paper. 

4.    Continuants  Applied  to  A  General  Ladder  Network. 

Let  Fig.  1  represent  a  general  ladder  network  constructed  of 

admittances  YOJ  Y.  ,  Y,  ,  .  .  .  Y  and  impedances  Z..  ,  Z„  ,  Z,. ,  ...  Z   ... 
2'   4'   6'     n  1    3'   5       n-1 
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Fig.  1.   General  Ladder  Network 

The  impedance,  Z  ,  of  the  network  looking  into  terminals  T„T, 
r         n  0  1 

can  be  written  as  a  continued  fraction  by  use  of  the  rules  for 
combining  series  impedances  and  shunt  admittances.   This  impedance  is 
given  by 


Z    -  Z,   + 
n     1 


Y2  + 


z3  + 


Y4  +   ' 


.   + 


Y   _   +   1 
n-2 


n-1  +  1 

Y 


n 


(4.1) 


which  can  be  expressed  as  a  ratio  of  continuants  as 
Z 


k(ZlfT) 


k(Y2,Yn) 


(4.2) 


Input  current,  I.  ,  to  the  ladder  network  is  given  by 


V 

z 


=  V 


k<Y2»  V 

k(Z  ,  Yn) 


(4.3) 
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Current  I-  is  that  fraction  of  I.,  given  by 


I3-  S 


Admittance  of  network  (Z0 , Y. }ZC , Yr  . . . .Z   n,Y  ) 

3'  4  5   6     n-1'  n 


Admittance  of  network  (Yn,Z,.Y.,Zc,,..Z   ,  .Y  ) 
,  2  3  4  5     n-1   n 


(4.4a) 


=  I. 


k(Z3>Yn) 


k(YvYn) 


2  k(Z3,Yn) 


Y„+ 


=   I. 


k<W 


Y2  k(Z3JYn)   +  k(Y4,Yn) 


(4.4b) 


Using  theorem  2  of  Appendix  I,  this  further  reduces  to 

k(vV 


i3  =  1 


In  a  similar  fashion 


MY2,Yn) 


(4.4c) 


I5     =     I3 


k<Y6>V 


k(Y4'V 


-     I. 


p     I. 


k(Y4,Yn)  k(Y6,Yn) 


Lk(Y2,Yn)  k(Y4,Yn) 


k<VV 


kTvTT 


(4.5a) 


(4.5b) 


In  general,  the  current  I»   _  through  any  series  impedance  Z_   1  can 


be  related  to  the  input  current  I,  as 


I2r+1~  *1 


k<Y2r+2'V 


k(Y2JYn) 


(4.6) 


by  a  continuation  of  the  above  process.   The  ratio  of  current  I 

through  Z„  . .  to  current  I~  , 1  through  Z„  . n  is  given  by 
2r+l  2m+l      °    2m+l    °  J 


2r+l 
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"214-1 

"2m+l 


T   k(Y0  ,0,Y  ) 
I     2r+2*  n^ 

k(Y2SYn) 
I  k(Y2m+2,V 


k(Y2,Yn) 


k<Y2r+2'V 


k<Y2m+2'V 


(4.7) 


The  voltages  across  the  shunt  admittances  should  be  considered 


next.   The  voltage  V     is  given  by 


V0   -  V 
2      g 


Impedance  of  network  (Y0,Z0 , Y, ,ZC , . . . .  Z   .  Y  ) 
r  2'  3'  4   5      n-1,  n 


(4.8a) 


=  V 


k(Z3,Yn) 


k(zryn) 


k(Y2,Yn)    •     k(Y2,Yn) 


(4.8b) 


=  V 


k(Z3?Yn)' 


(4.8c) 


In  similar  fashion  the  voltage  V,  is  given  by 

fk<Z3>V 


\     •     V2 


k(Z5'V 
k<Z3'V 


k<zi-V 


k(Z5'V 
k<Z3>V 


(4.9a) 


-  V 


k(Z5'V 
k(Z  YH) 


(4.9b) 


In  general,  the  voltage  V„   across  shunt  admittance  Y0   can  be 
&      j  &   2r  2r 


related  to  the  input  voltage  V  as 

g 


V„   =  V 
2r      g 


k(Z0  ,.  ,Y  ) 
2r+l   n7 


k(Z1?Yn) 


(4.10) 


The  ratio  of  voltage  V0   across  shunt  admittance  Y„   to  voltage  V0 

°   2r  2r         °        2m 

across  shunt  admittance  Y„   is  given  by 
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V 


2r 
2m 


k(Z0  .. ,Y  ) 
2r+l-  n 


k(Z0  .,  ,Y  ) 
2m+l   n  „ 


(4.11) 


It  is  often  of  interest  to  determine  the  voltage-ratio  transfer 

function,  V  /V  ,  of  the  ladder  network.   This  can  be  done  in  the 
n  g' 

following  manner: 
letting 

2r+l  =  n-1 
then 

n  -  2r+2 
Using  Eq.  (4.6),  we  get 


■n-1 


=  I. 


k(Y  ,Y  ) 
n  n 


k(Y2,Yn) 


(4.12) 


Since 


V  /Z 
g  n 


(4.13) 


then 


V 


n-1 


k(Y  ,Y  ) 
n  n 


Y  k(Y„,Y  ) 
n   2'  n 


V 

Z 

n 


(4.14) 


and  substituting  from  Eq.  (4.2)  for  Z  gives 


V 


k(Y   ,Y  ) 
n     n 


Y  k(Y0,Y   ) 
n        2      n 


k<Y2-V 
k(Zi;Yn) 


[ 


k(Yn'V 
Yak<Zl'V 


(4.15) 


but 


k(Y    ,Y   )      =     Yn 
n     n  n 
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therefore 


V 
_r 

v" 


k(Z15Yn) 


(4.16) 


If  the  ladder  network  is  such  that  Z1  =  0,  then  using  Theorem  6  of 


Appendix  I  gives 


k(Z3,Yn) 


'n     k(Y2JYri) 


(4.17) 


Vg     k(Z35Yn) 


(4.18) 


and 


V 


2r 


=  V 


=  V 


k<Z2r+l'V 


■2r+l   = 


k(Z3,Yn) 


k<W 
k(Z3,Yn) 


k(Y2r+2,Yn) 


k(Y4'V 


(4.19) 


(4.20) 


(4.21) 


5.   First  and  Second  Order  RC  Ladder  Networks. 

It  is  of  interest  to  a  design  engineer  to  be  able  to  synthesize  an 
electric  network  from  a  given  voltage-ratio  transfer  function.   The 
problem  is  that  of  determining  the  values  of  circuit  components  such 
that  they  will  provide  the  desired  response  when  connected  together  to 
form  a  circuit.   A  method  for  determining  these  values  for  a  first  Cauer 
form  of  RC  ladder  network,  such  as  shown  in  Fig.  2.,  will  now  be  illus- 
trated by  starting  with  a  single  resistor-capacitor  network  and  then  . 
progressively  adding  on  more  components. 
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Fig.  2.   RC  Ladder  Network  of  the  First  Cauer  Form 
The  voltage-ratio  transfer  function  of  the  network  as  given  by  Eq. 
(4.16)  as  a  function  of  the  complex  frequency  variable,  s,  is 


V 


Q  (s)     k(R. ,sC  ) 
m  1   n 


(5.1) 


For  a  two-component  network  the  denominator  polynomial  is 

Q2(s)      =      R^s+1      =      Tns+1 
Define   a  L-section   time   constant,    T. . ,    as 


T..   =  R.C. 
li      li 


(5.2) 


(5.3) 


for  various  values  of  the  subscripted  variables.   If  we  are  given  a 
transfer  function  whose  denominator  polynomial  is  of  the  form 

(5.4) 
where  K,  is  some  positive  real  constant,  and  we  wish  to  synthesize  a 


Q2(s)   =   KLs+l 


network  for  this  function,  then  we  can  set 
RnCn   =  K, 


and  get 


11 


R, 


Kn 


(5.5) 


By  assigning  a  positive  value  to  C  ,  we  can  calculate  a  value  for  R  , 
and  hence  will  have  the  necessary  component  values  to  synthesize  the 
network. 
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Thus  far  this  has  been  a  rather  trivial  problem.   Adding  on  an- 
other resistor  and  capacitor  will  complicate  the  situation  somewhat. 

For  a  four -component  network  the  denominator  polynomial  of  the 
transfer  function  is 

Q4(s)   =  TuT22s2  +   (Tu  +  T12  +  T22)s  +  1  (5.6) 

The  denominator  polynomial  of  a  given  transfer  function 
would  be  of  the  form 

Q4(s)   =  KLs2  +  K2S  +  1  (5.7) 

where 

Kl  =  T11T22 

K2  =  Tll  +  T12  +-  T22  (5*8) 

From  Eqs.  (5.8)  we  get 

T12  =  K2   '   Tll   "  T22  (5,9) 

Let 

Tn   =  a  =  R1C1 

where  a  is  some  positive  real  constant.   (Note:   negative  values  of  a_ 
have  no  meaning  since  this  would  imply  a  negative  time  constant  in  the 

ladder  network).   Eq.  (5.9)  can  be  rewritten  as 

Kl 
T12  =  K2  "  a  *  a"  ,=  R1C2  (5*10) 

Define  T. .  as  a  cross  time  constant.   It  is  the  product  of  R.  and  C. 

for  all  values  of  the  subscripted  variables  except  for  i  equal  to  j. 

All  time  constants  for  the  network  must  be  positive  since  we  consider 

only  positive  values  for  the  R's  and  C's.   The  requirement  that  T1 «>  0 

will  be  satisfied  if  the  values  of  a  are  restricted  to  those  positive 

values  for  which 

Kl 
K0  >  a  +  —  (5.11) 

2        a  v 
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We   can  define   a   function,    f(a),    from  Ineq.    (5.11)    as 

Kl 
f(a)      =     K2   -   a   -  -±     >  0  (5.12) 

A  typical  plot  of  the  function  is  shown  in  Fig.  3.  for  assumed  values 
of  K,  equal  to  ten  and  K„  equal  to  twenty.   Since  we  have  stated  that 
f(a)  must  be  greater  than  zero,  then  the  graph  can  be  used  to  determine 
those  values  of  a  for  which  the  restriction  on  f(a)  is  satisfied. 

Using  Eqs.  (5.8)  and  (5.10)  and  the  time  constant  definition  of 
T1 , ,  we  can  write  the  following: 


C2 

= 

K2   - 

a   - 

a 

Ri 

Cl 

= 

a 

Ri 

R2 

= 

aC2 

(5.13) 


By  assigning  a  positive  value  to  R..  and  using  Eqs.  (5.13),  the  values 
of  C1 ,  C„,  and  R~  can  be  calculated  to  synthesize  the  desired  network. 
To  illustrate  with  an  example,  assume  that  we  are  given  a  voltage- 
ratio  transfer  function  whose  denominator  polynomial  is 

Q,(s)   =   10s2  +  20s  +   1 
where 

K-   =   10 

K2  =   20 
Fig.  3.  was  plotted  for  these  values  of  K,  and  K„,  therefore  it  can  be 
used  to  determine  an  allowable  value  for  the  constant,  a_.   We  see  from 
Fig.  3,  that  all  values  of  a  from  about  0.5  to  19.5  result  in  f(a)  being 
greater  than  zero,  hence  these  are  all  allowable  values.   If  we  choose 
a  to  be  unity  and  assign  a  value  of  9  megiohms  to  R1  ,  then  Eqs.  (5.13) 
will  give 
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Fig.  3.   Range  of  Allowable  Values  of  Time  Constant  a. 
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20-1-10       -6  .    , 
=   10    farads 


9  x  10 
R_   =  r—  =   10  megohms 

2    io-6 

C,   =  t     =  77  x  10    farads 

1     9  x  106     9 

The  component  values  can  be  substituted  into  Eqs.  (5.8)  to  insure  that 
they  do  in  fact  give  the  desired  transfer  function. 
6.   Third  Order  RC  Ladder  Networks. 

Synthesizing  a  two-component  ladder  network  was  seen  to  be  a 
trivial  problem.   Adding  two  more  components  to  the  network  made  the 
problem  a  little  more  complex,  but  still  quite  straight-forward. 
Adding  a  third  resistor  and  capacitor  to  form  a  six-component  network 
will  now  be  discussed  at  length  to  better  illustrate  the  synthesis 
procedure. 

For  a  six-component  network  the  denominator  polynomial  of  a  given 
voltage-ratio  transfer  function  is 

Q6(s)   =  KLs   +  K  s2  +  K3s  +  1  (6.1) 


where 


Kl     T11T22T33 


K2  =  T11T22  +  T11T23  +  T11T33  +  T12T33  +  T22T33     (6,2) 


K3     Tll  +  T12  +  T13  +  T22  +  T23  +  T33 


Let 


Tll  =  a  and  T33  =  b 
where  a  and  b  are  some  positive  real  constants.   Eqs.  (6.2)  can  be 
rewritten  as 


25 


Kx   -  ab  T22  (6.3a) 

K2  =  aT22  +  aT23  +  ab  +  bTi2  +  bT22           (6.3b) 

K3  =  a  +  T  +  T13  +  T22  +  T23  +  b            (6.3c) 
Solving  for 

Ki 

T22   "  Sb  (6-M) 

from  the  first  of  Eqs.  (6.3)  and  substituting  this  into  the  second 
two  gives 


bTu 

+ 

aT23 

:  K2 

-  ab   • 

Ki 

b 

a 

(6.4a) 

T12 

+ 

T 
13 

+ 

T23 

-  K3 

a 

-  b 

-  Kl 
ab 

(6.4b) 

If  the  network  is  to  be  physically  realizable,  then 

K      K 

K_   -   ab   -  —-  -  —  >  0  (6.5a) 

2.  b      a 

Kl 
KQ   -   a   -  b   -  ~   >0  (6.5b) 

3  ab 

must  hold  since  all  values  of  the  R's  and  C's  are  to  be  positive.   The 
left  side  of  Eqs.  (6-4)  are  therefore  greater  than  zero.   Ineqs.  (6.5) 
can  be  considered  as  necessary  conditions  on  allowable  values  for  a 
and  b_.   It  should  be  noted  that  these  inequalities  will  hold  if  the 
given  transfer  function  was  derived  from  an  actual  RC  ladder  network 
of  the  first  Cauer  form.   (This  can  be  proven  by  using  simple  in- 
equalities and  the  fact  that  we  started  out  with  all  positive  values 
in  Eqs.  (6.3)). 

Investigation  of  Ineqs.  (6.5)  should  prove  to  be  instructive  at 
this  point  since  they  were  stated  to  be  necessary  conditions  for 
physical  realizability  of  the  network.   If  we  set  these  inequalities 
equal  to  zero  and  solve  for  a  in  terms  of  b  for  both  equations,  then 

we  can  plot  a   as  a  function  of  b  for  given  values  of  K..  ,  K9,  and  K«. 
The  plots  will  show  allowable  values  for  a  and  b. 
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Setting  Ineqs.  (6.5)  equal  to  zero  and  solving  for  a  in  terms  of 

b  gives  , 

(I^-I^b)   +   Jo^b-K^)2  -  4Kxb3 

a  =  t-x (6.6a) 

-2,b 

from  Ineq.  (6.5a)  and 


(b  -K~b)   +   J(K„b-b2)2  -4K,b 

a  =  - —  (6.6b) 

-2b 

from  Ineq.  (6.5b).   If  we  had  solved  for  b_  in  terms  of  a_,  we  would 

have  gotten 


4  -   K2a)  ±  \| 
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(K  -  K  a)   +  J  (K„a  -  K  )   -  4K  a 

Y (6.7a) 

-2a 


from  Ineq.  (6.5a)  and 

2  1         9  2 

(a  -  K„a)  +   J  (K  a  -  a  )      -  4K, a 

b  =  ^ ^ — L—  (6.7b) 

-2a 

from  Ineq.  (6.5b).   Note  the  symmetry  between  Eqs.  (6.6)  and  (6.7). 

Allowable  values  of  a  and  b  must  be  both  positive  and  real,  hence  the 

plots  are  in  the  first  quadrant.   The  allowable  values  are  restricted 

to  those  values  which  lie  within  the  region  bounded  by  both  of  the 

closed  curves.   Plots  of  Eqs.  (6.6a)  and  (6.7a)  result  in  identical 

curves  because  of  the  symmetry  of  the  equations.   The  same  is  true  of 

plots  of  Eqs.  (6.6b)  and  (6.7b). 

Fig.  4.  shows  the  curves  of  Eqs.  (6.6)  for  an  assumed  given  case  of 

Q6(s)   =  3s3  +  9.5s2  +  8s  +  1  (6.8) 

where 

KL   =  3      K2  =  9'5      K3   =  8 
Eqs.  (6.4)  serve  as  a  basis  in  attempting  to  synthesize  the  network. 
Suitable  values  for  the  constants  a   and  b  must  be  chosen,  or  be  given. 

Three  interesting  cases  arise  in  considering  the  problem.   These  cases 
are: 
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Fig.  4.   Range  of  Allowable  Values  of  Constants  a  and  b 
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A)  We  are  free  to  choose  values  of  a  and  b  so  long  as  Ineqs.  (6.5) 
are  satisfied. 

B)  We  are  given  a  specific  value  of  a  or  b  and  must  choose  the 
other  such  that  Ineqs.  (6.5)  are  satisfied. 

C)  We  are  given  specific  values  for  both  a  and  b. 
Each  of  these  three  cases  will  be  discussed  separately. 

A.   Constants  a  and  b  Not  Specified 
If  a  and  b  are  not  specified,  then  we  are  free  to  choose  them  to  be 
any  value  desired  so  long  as  Ineqs.  (6.5)  are  satisfied  and  the  final 
circuit  component  values  satisfy  Eqs.  (6.2).   At  this  point  it  is  con- 
venient to  let  a  equal  b  in  order  to  simplify  subsequent  equations. 
(Justification  for  assuming  that  a  equals  b  lies  in  the  symmetry  of 
Eqs.  (6.6)  and  (6.7).   If  some  value  of  a  is  allowable  then  certainly 
that  same  value  of  b  is  also  allowable  as  far  as  Ineqs.  (6.5)  are  con- 
cerned.)  Eqs.  (6.4)  can  then  be  rewritten  as 

K2  2KX 

T12  +  T23   =  T  '      B      '     ~1  (6'9a) 

a 

Kl 
T12  +  T  3  +  T23   =  K3  -  2a   -  -  (6.9b) 

a 
Subtracting  the  first  equation  from  the  second  gives 

T13   -  K3  -J-     -  a  +  ^  (6.10) 

a 

Since  T, ^  must  be  a  positive  real  constant,  then 

K        K 

K3  "  r  "  a  +  ~2   >0  <6ai) 

a 

must  also  be  satisfied.   Ineq.  (6.11)  can  be  considered  as  a  third 

necessary  condition  on  the  allowable  values  of  a.   Those  values  which 
satisfy  Ineqs.  (6.5)  must  also  satisfy  Ineq.  (6.11)  in  order  to  be 
allowable.   A  plot  of  Ineq.  (6.11)  as  a  function  of  a  for  given  values 


29 


of  K1 ,  Kn,  and  K-  can  be  made  to  find  the  values  which  are  allowable. 
(As  an  alternative,  values  of  a_  which  will  satisfy  Ineq.  (6.11)  are 
found  on  the  segment  of  a  line  of  unit  slope  through  the  origin  of  the 
plots  of  Eqs.  (6.6)  that  is  bounded  by  the  innermost  curves).   Any 
allowable  value  of  a  can  be  used  in  the  synthesis  procedure. 

After  determining  a  suitable  value  of  a_,  we  can  calculate  T..  „  from 
Eq.  (6.10)  for  the  given  values  of  K..  ,  K„,  and  Ko.   By  definition, 

T13   =  R1C3 

or 

T13 
C3   -  ^  (6.12) 

Assigning  a  positive  value  to  FL  allows  us  to  calculate  CL.   C,  and  R„ 

can  be  found  from 

C   =  — 
1     R, 


and 


R3     C3 


since  a  =  b  and  by  definition 


Tll  =  R1C1   =  a 


and 


T33   ~  R3C3   "  b 
Note  that  in  effect,  Eq.  (6.12)  gives  a  relationship  between  time  con- 
stants T, ,  and  Too  by  giving  a  relationship  between  R..  and  Co. 

We  have  thus  far  found  suitable  values  for  four  of  the  six  circuit 
components.   We  must  still  find  values  for  R_  and  C„  in  order  to  synthe- 
size the  network.   These  can  be  determined  from  Eq.  (6.9a)  and  the  fact 

that  K. 

22     *2^2      2 
a 
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In  Eq.  (6.3d).   Substituting 

Kl 
C   =  -±-r  (6.13) 

R2aZ 

into  Eq.  (6.9a)  and  using  time-constant  definitions  give 
R,  K,  K~  i 

fi  +  R2C3  -  -t  "  a  "  -r  (6-14) 

R2a  a 

from  which  we  can  obtain  an  expression  for  R„  by  using  the  Quadratic 
Formula.   Then 

"~1 


*/ 


4K  T 
Y2 


R2  -  ^ *- (6.15) 

where  Y  is  defined  as 

K9  2K1 

Y  =  a   "  a  "    2  (6,16) 

a 

as  a  convenience  in  writing  Eq.  (6.15).   Knowing  values  for  a,  K. ,  K„, 
R,  ,  and  CL  allows  us  to  solve  for  R  .   Finally,  we  can  find  C„  from  Eq. 
(6.13)  and  our  work  is  finished. 

In  calculating  R  using  Eq.  (6.15),  it  should  be  noted  that  only 
those  values  of  a  which  were  previously  determined  to  be  allowable  can 
actually  be  used  in  the  synthesis  procedure  for  which 

a2Y2  -  4K1R1C3^  0  (6.17) 

Since  we  cannot  have  a  complex  resistance.   Ineq.  (6.17)  represents  a 
fourth  restriction  on  allowable  values  of  a.   It  can  be  plotted  as  a 
function  of  a  to  determine  which  previously  allowable  values  are  still 
allowable. 

We  have  in  effect  then  been  able  to  take  Eqs.  (6.2),  which  are  three 
equations  in  six  unknowns,  and  solve  for  the  six  unknowns.   This  was  done 
by  finding  suitable  values  for  two  time  constants  a   and  b  (where  a  was 
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assumed  to  be  equal  to  b  )  ,  assigning  a  value  to  one  of  the  unknowns,  R.  , 
and  solving  for  the  other  five  unknowns  in  terms  of  R  and  a  such  that 
Eqs.  (6.2)  are  satisfied. 

It  might  seem  at  first  glance  that  finding  values  of  a  which 
satisfy  the  four  restrictions  given  by  Ineqs.  (6.5),  (6.11),  and  (6.17) 
would  prove  to  be  so  laborious  as  to  make  the  synthesis  procedure  of 
little  value.   Such  is  not  the  case,  especially  if  one  has  access  to  a 
digital  computer.   (Making  the  necessary  plots  by  hand  might  require  a 
bit  of  work,  but  the  task  is  quite  straight-forward.)   It  should  be 
noted  from  Eqs.  (6.15)  and  (6.13)  that  for  each  allowable  value  of  a 
and  each  assigned  value  of  R  ,  we  get  two  sets  of  values  of  R„  and  CL 
which  can  be  used  to  synthesize  the  network.   It  should  be  noted  that 
there  generally  will  be  any  number  of  allowable  values  of  a   for  some 
given  transfer  function,  hence  any  number  of  different ■ sets  of  component 
values  that  can  be  found  to  do  the  job. 

As  an  example  of  the  procedure,  assume  that  we  have  a  given  trans- 
fer function  whose  denominator  polynomial  is  the  same  as  that  given  by 
Eq.  (6.8).   The  plots  of  Ineqs.  (6.5)  for 
Kl  =  3     K  =  9.5     K3  =  8 

shown  in  Fig.  4.  can  be  used  in  determining  allowable  values  for  a_.   A 

convenient  choice  for  a   might  be  unity  since  this  would  help  reduce  the 

mathematical  calculations,  and  unity  is  an  allowable  value  as  shown  by 

Fig.  4.   Using  Eq.  (6.10),  we  have 
K  K 

T13  =  K3  "  ~T~  '      l   +  1    =   (8-9-5"1+3)  =  °-5 
If  we  let 

R   =  0.5  ohm 
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then 

T 


R^  =  S3  =  l    farad 


By  using  the  initial  assumption  that 


a   =   Tn   =  RlCl 


and 


b  "  T33   "  R3C3   "  a 


we  can  determine  that 


ci  =  t:  =  h  =  2  farads 


and 


a      1     n   , 
R~   =  —  =  -j-  =   1  ohm 

5  S     L 


Using  Eq.  (6.16),  we  find  that 

Y  -   (9.5  -  1  -  6)   =  2.5 
Eq.  (6.15)  gives 


2.5  +  V  (2.5)2  -  4(0.5)3  =  2.5  ±  0.5        _ 

R.   =  = v  ■  - —  = =   1;  1.5  ohms 

2  2  2 

For  R~  equal  to  one  ohm, 

Kl      3 
C„   =  j  =  y  =  3  farads 

R2a 
and  for  R_  equal  to  one  and  one-half  ohms, 

3 
Cy     =     - — -     =     2   farads 

The  component  values  for  the  synthesized  network  are  then: 
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Network  #1  Network  #2 

R   =  0.5  ohm  ^   =  0.5  ohm 

R„   =   1    ohm  R„   =   1.5  ohms 

R„   =   1    ohm  ,  R„   =   1    ohm 

3  and  3 

C   =  2   farads  C   =   2    farads 

C       =  3    farads  C       =  2   farads 

C3   =   1   farad  C,   =   1    farad 

As  a  check  for  Network  #1: 

Kl   =   T11T22T33   =   3 

K=TT  +TT  +TT  +TT  +TT 

K2  Hi    22  Hi   23  illi33  i12i33  i22i33 

=      3      +      1      +      1     +      1.5     +     3      =      9.5 

K3  =  Tll  +  T12  +  T13  +  T22  +  T23  +  T33 
=   1  +  1.5  +  0.5  +  3  +  1  +  1  =  8 
As  a  check  for  Network  #2: 
Kt  =  3 

K2  =  3  +   1. 5  +  1  +  1  +  3   =  9. 5 

K3   =   1  +  1  +  0.5  +  3  +   1.5  +   1   =  8 

These  values  bf  K  ,  K„,  and  K-  check  with  those  of  the  given  transfer 
function,  and  Eqs.  (6.2)  are  satisfied. 

B.   Constants  a  or  b  Specified 
A  situation  might  arise  in  which  one  wishes  to  synthesize  a  net- 
work from  a  given  voltage-ratio  transfer  function  with  a  particular 
value  of  a  or  b  to  be  used.   This  is  equivalent  to  saying  that  we  wish 
the  network  to  have  a  particular  value  for  the  time  constant  of  the 
first  or  last  L-section  of  the  network.   The  logical  thing  to  do  would 
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seem  to  be  to  let  a  equal  b  again  and  see  if  this  value  for  a  satisfies 
Ineqs.  (6.5),  (6.11),  and  (6.17).   If  so,  then  the  procedure  outlined 
in  part  A  above  can  be  used  to  synthesize  the  network.   If  not,  then 
the  procedure  outlined  in  part  C  to  follow  can  be  used. 
C.   Constants  a_  and  b  Specified 
If  we  are  required  to  use  particular  values  of  a   and/or  b,  or  if 
we  did  not  wish  to  assume  that  a  equals  b_,  then  a  more  general  ap- 
proach must  be  used  to  synthesize  the  network.   We  must  fall  back  on 
Eqs.  (6.4)  as  a  basis  for  developing  a  procedure.   If  the  specified 
values  of  a  and  b_  do  not  satisfy  Ineqs.  (6.5),  we  cannot  physically 
realize  the  network.   In  addition,  if  either  a   or  b  is  specified  and  we 
cannot  find  a  value  for  the  other  such  that  Ineqs.  (6.5)  are  satisfied, 
then  again  we  cannot  physically  realize  the  network. 

Assume,  however,  that  the  values  of  a  and  b_,  either  specified  or 
chosen,  do  satisfy  Ineqs.  (6.5).   The  problem  now  is  to  synthesize  the 
network  using  these  values. 

It  is  convenient  at  this  point  to  define 

K      K 
W  =  K2  -  ab  -  -  -  - 

and 

Ki 

Z  "  K3  "  a  "  b  "  3b 

where  W  and  Z  are  positive  real  constants  if  Ineqs.  (6.5)  are  satisfied. 
Eqs.  (6.4)  can  be  rewritten  in  more  compact  form  as 

bT12  +  aT23  =  W  (6.18a) 


T12  +  T13  +  T23   =  Z  (6.18b) 
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The  problem  now  is  to  first  solve  Eqs.  (6.3a)  and  (6,18)  for 
suitable  values  of  R  ,  C        R  ,  and  CL .   The  procedure  will  be  to  solve 
for  three  unknowns  in  terms  of  the  fourth,  and  then  assign  a  suitable 
value  to  that  fourth  unknown.   After  determining  values  for  R  and  C- , 
we  can  readily  obtain  values  for  C.  and  R~  from  the  time  constant 


definitions 


and 


Tll   "   R1C1   "   a 


T33   "  R3C3   "  b 
(It  should  be  noted  that  C1  and  R„  do  not  appear  in  the  basic  Eqs. 
(6.3a)  and  (6.18).   The  only  constraint  imposed  on  the  values  of  C.  and 
R~  is  that  the  product  of  R1  and  C.  must  equal  the  allowable  value  of  a 
and  the  product  of  R„  and  CL  must    equal  the  allowable  value  of  b). 
Eqs.  (6.3a)  and  (6.18)  can  be  rewritten  as 

C2  "  Sq  <6'19a> 

a  W 

RnCL  +  r  RoC-3  =  -  (6.19b) 

1  I  b   z  J     a 

R  Cl     +     R,C3  +  R2C3   =  Z  (6.19c) 

Substituting  the  expression  for  CL  from  the  first  equation  into  the  last 
two  gives 

R  K 

ihq  +  r-R2c3  =  !  <6=20a> 

RiKi 

-^       +     R1C3  +  R2C3   =   Z  (6.20b) 

Eqs.  (6.20)  are  two  equations  in  three  unknowns.   Solving  Eq.  (6. 20a) for 
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a 

-  RiKi 

abR2 

—  R 
b  R2 

(6.21) 


and  substituting  into  Eq.  (6.20b)  gives 


abR„ 


+ 


_W 
a 


abR„ 


^  R 
b  R2 


+  R, 


a 

.  R1K1 
abR2 

a 
b 

-r2 

=  z 


(6.22) 


which  is  one  equation  in  two  unknowns.   Applying  the  Quadratic  Formula 
to  Eq.  (6.22)  results  in 


R2  -  Rx 


(abK   +  Wb   -   b2K  ) 

2~2        3 
2(a  b  Z  -  Wb  ) 


(6.23) 


+  ^  (abKL  +  Wb3   -   b2Kx)2  +  ^^(Wb3  -  a2b2Z) 


2  2        3 
2 (a  b  Z   -  Wb  ) 


The  quantity  in  brackets  in  Eq.  (6.23)  can  be  calculated  from  the  known 
constants  a,  b,  K  ,  K?,  and  K-.   It  should  be  noted  however  that  only 
the  allowable  values  of  a  and  b_  derived  from  the  plots  of  Ineqs.  (6.5) 
which  result  in  the  quantity  under  the  radical  sign  being  equal  to  or 
greater  than  zero  can  be  used.   Assigning  a  suitable  value  to  R1  allows 
us  to  calculate  R„.   Then  determine  C„  from  Eq.  (6.21).  The  last  three 
component  values  can  be  obtained  from 

c    =  -4- 


R, 


and 


b 
C3 

A 
abR, 
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This  then  gives  values  for  the  six  components  of  the  network.   Note 
that  again  we  get  two  values  for  R~  from  Eq.  (6.23)  and  hence  two 
values  of  C-  from  Eq,  (6. 19a)  for  each  set  of  allowable  values  of  a 
and  b_  for  every  assigned  value  of  R-  . 

A  slightly  different  approach  was  used  for  finding  R„  in  the  last 
case  as  compared  to  that  used  for  the  case  where  a  =  b.   However,  the 
results  are  the  same,   This  can  be  seen  by  substituting  a  =  b  into  Eq. 
(6.23)  and  using  the  definitions  of  W,  Y,  and  Z  and  Eq .  (6.12),   The 
approach  used  in  the  last  case  is  more  general,  as  we  shall  subsequently 
see, 
7.    Fourth  Order  RC  Ladder  Networks. 

The  ideas  expressed  and  techniques  developed  in  synthesizing  first, 
second,  and  third  order  RC  ladder  networks  can  be  extended  to  higher 
order  systems.   This  can  be  illustrated  by  adding  another  RC  L-section 
to  form  a  fourth  order  network.   The  first  thing  noticed  is  that  the 
algebraic  expressions  are  getting  longer  and  more  involved.   Neverthe- 
less, we  can  proceed  in  the  same  manner  as  before. 

For  an  eight-component  network  the  denominator  polynomial  of  a 

given  voltage-ratio  transfer  function  is 

4        3        2 
Q8(s)   =  K  s   +  K2s   +  K3s   +  K.s  +   1  (7.1) 

where 

Kl  =  T11T22T33T44  (7,2a) 

V  =TfTpr  imrprp  4-TTT  4-rTTirp 

2     illi22i33     illi22134     illi22i44    11  2,3  44 


+  T11T33T44  +  T12T33T44  +  T22T33T44  (7'2b) 
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K3  "  T11T22  +  T11T23  +  T11T24  +  T11T33  +  T11T34 


+  T11T44  +  T12T33  +  T12T34  +  T12T44  +  T13T44 


+  T22T33  +  T22T34  +  T22T44  +  T23T44  +  T33T44 


(7.2c) 


K4  "  Tll  +  T12  +  T13  +  T14  +  T22  +  T23  +  T24 

+   T     +   T     +   T 

^3     x34      44  (7. 2d) 


Let 


Tll  =  a       T44  =  a        T22  =  b 

Setting  T. .  and  T,,  equal  to  the  same  constant  is  justified  in  that  the 

only  restriction  on  values  of  Cn  and  R.  is  that  the  products  R, C,  and 

14  11 

R,C,  be  equal  to  a.   The  components  C  and  R,  do  not  show  up  in  any  of 
the  cross  time  constants,  T...   After  finding  proper  values  for  R..  and 
C,  ,  we  can  always  adjust  C,  and  R,  so  that  their  respective  products  do 
equal  a. 

With  the  L-section  time  constants  thus  defined,  we  can  rewrite  Eqs, 

(7.2)    as 

Kl 
T33     =     —  (7.3a) 

a  b 
aiT12  +aS3  +     abT34  =     X  (7'3b) 


aLb 


(-\-  +a)T12+aT13+2aT23+aT24+(a+b)T34+T12T34  =  Y  (7.3c) 

T12  +  T13  +  T14  +  T23  +  T24  +  T34   =  Z    (7'3d) 


where 


2K  K 

X  =  K0 -  a  b  -  —■  >  0 

2   a  b 
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Y  =  K3  -  2ab 


2    2K1 
1   "  aT" 


K, 


K 


1 


■i  *  ° 

a 


(7.4) 


Z  =  K,  -  2a  -  b 7  >  0 

4  2U 

a.  b 
The  requirement  that  X,  Y,  and  Z  all  be  greater  than  zero  must  be 
satisfied  if  the  network  is  to  be  physically  realizable.   Ineqs. 
(7.4)  can  be  considered  as  necessary  conditions  on  the  allowable  val- 
ues of  a   and  b  for  a  given  transfer  function.   Again  plots  of  Ineqs. 
(7.4)  can  be  made  to  determine  allowable  values  of  a   and  b_  for  use  in 
the  synthesis  procedure. 

Initially  we  started  out  with  Eqs.  (7.2),  which  are  four  equations 
in  eight  unknowns.   By  finding  suitable  values  for  a  and  b_,  we  are  able 
to  write  Eqs.  (7.3),  which  are  four  equations  in  six  unknowns.   The 
number  of  unknowns  in  the  equations  can  be  further  reduced  by  noting  the 

following  cross  time  constant  definitions: 

b 


-L-i  ry  K-i  yjm  i\... 


T13  "  R1C3  ~  Rl 


T23     R2C3     R2 


R, 


In  the  above,  C„  and  C,  have  been  expressed  in  terms  of  R„ ,  R„ ,  and 
known  constants.   The  last  three  of  Eqs.  (7.3)  can  be  rewritten  as 
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K, 


R, 


K, 


K 


R    + 


1 
ab 


m 


+   abR0C, 
3  4 


R, 


2K  /R  \ 

ab-(R;)+aR2V<a+b)R3C4 


2  4     J  4     14 


=  X    (7.5a) 


=   Y    (7.5b) 


=   Z 


(7.5c) 


which  are  three  equations  in  four  unknowns.   Then  proceed  to  solve  for 
three  of  these  unknowns  in  terms  of  the  fourth.   It  should  be  noted 
that  only  one  component  of  each  of  the  four  L-sections  of  the  network 
shows  up  in  Eqs.  (7.5)  as  one  of  the  unknowns. 

Solving  Eq.  (7.5a)  for  C,  and  substituting  it  into  Eqns.  (7.5b) 


and    (7.5c)    gives 


—  A 

a         3 


•2  t  [< 


i2b2+bX-K1)R2R3    "   K1R2  R3 


R1  + 


4 
KxaR2 


K  a 


+      (Kj+aX   •     ~  )R2R3     +      (aX+bX-abY)R2R3  =        ° 

Kl  2  Kl    2  Kl        2 

"7*3  Rl     +      ["  ~R2     +  XR2R3     +      (ab-   a"  )R3   J     Rl 

R2R2J      =     0 


(7.6a) 


\3 
b   R2 


1  9 

(—  +X-  ^RoRq     +      (X-abZ): 
v  a-  b        2  3  v 


(7.6b) 


which  are  two  equations  in  three  unknowns.   Eq.  (7.6a)  is  quadratic  in 
R   .   It  can  be  solved  for  R  by  use  of  the  Quadratic  Formula.   This 
will  give  two  values  of  R.  to  be  put  into  Eq.  (7.6b),  one  at  a  time,  to 
obtain  one  equation  in  the  two  unknowns  R_  and  R».   Admittedly  this  last 
equation  will  be  rather  complicated.   However,  we  are  able  to  assign  a 
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value  to  either  R  or  R„  and  solve  for  the  other,  graphically  if  neces- 
sary since  this  is  equivalent  to  changing  the  impedance  level  of  the 
network.   After  finding  values  for  R  and  R„ ,  we  can  go  back  and  find 
corresponding  values  for  the  other  six  components  of  the  network. 
8.    Nth  Order  RC  Ladder  Networks. 

By  now  it  should  be  evident  that  a  procedural  pattern  has  been 
evolved  for  the  synthesis  procedure.   We  can  find  any  number  of  sets 
of  component  values  which  will  satisfy  a  given  voltage-ratio  transfer 
function,  but  these  component  values  will  be  restricted  to  certain 
ranges.   The  ranges  will  be  determined  by  the^range  of  allowable  values 
of  the  L-section  time  constants.   It  should  also  be  obvious  that  a  digi- 
tal computer  is  needed  if  we  are  to  synthesize  large  order  networks. 

An  analysis  of  what  has  taken  place  in  the  synthesis  of  RC  ladder 
networks  up  to  the  fourth  order  shows  more  clearly  the  procedural  pat- 
tern, and  how  the  procedure  might  be  extended  to  the  synthesis  of  an  nth 
order  system. 

In  the  first  order  network  there  was  an  initial  system  equation  in 
two  unknowns.   This  was  Eq.  (5.2).   There  were  no  necessary  conditions 
involved  and  no  allowable  values  of  L-section  time  constants  to  be  deter- 
mined.  By  assigning  a  value  to  one  unknown  in  Eq.  (5.5),  we  could  solve 
for  the  other. 

In  the  second  order  network  there  were  two  initial  system  equations 
in  four  unknowns.   These  were  Eqs.  (5.8).   We  had  one  necessary  condition 
to  satisfy  (Ineq.  (5.12))  and  one  L-section  time  constant  for  which  to 
determine  allowable  values.   With  this  time  constant  we  were  able  to 
reduce  the  two  initial  system  equations  to  one  equation  in  two  unknowns. 
Assigning  a  value  to  one  of  the  unknowns  allowed  us  to  find  values  for 
the  other  three  using  Eqs.  (5.13). 
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The  third  order  system  had  Eqs.  (6.2)  as  three  initial  system 
equations  in  six  unknowns.   Ineqs.  (6.5)  were  two  necessary  conditions 
which  had  to  be  satisfied  if  the  network  was  to  be  physically  realized. 
It  was  necessary  to  determine  allowable  values  for  two  L-section  time 
constants  which  would  satisfy  these  conditions.   It  was  also  noted  that 
the  two  time  constants  could  be  made  equal  to  each  other  if  desired. 
The  time  constants  were  used  to  reduce  the  three  initial  system  equa- 
tions to  Eqs.  (6.3d)  and  (6.4),  which  are  three  equations  in  four  un- 
knowns.  We  then  solved  for  three  of  these  unknowns  in  terms  of  the 
fourth,  assigned  a  value  to  the  fourth,  and  subsequently  went  back  and 
found  values  for  all  the  others. 

In  the  fourth  order  system,  Eqs.  (7.2)  were  four  initial  system 
equations  in  eight  unknowns.   There  were  three  necessary  conditions 
given  by  Ineqs.  (7.4)  to  satisfy.   It  was  necessary  to  determine  allow- 
able values  for  three  L-section  time  constants  which  would  satisfy  these 
inequalities.   It  was  noted  that  one  can  always  set  the  first  and  last 
L-section  time  constants  equal  to  each  other.   The  time  constants  enabled 
us  to  reduce  the  initial  system  equations  to  Eqs.  (7.3),  which  are  four 
equations  in  six  unknowns.   By  redefining  some  of  the  cross  time  con- 
stants, we  further  reduced  the  system  equations  to  Eqs.  (7.5)  and  Eq. 
(7.3a),  which  are  four  equations  in  five  unknowns.   We  again  solved  for 
four  of  these  unknowns  in  terms  of  the  fifth,  assigned  a  value  to  the 
fifth,  and  went  back  to  find  values  for  all  the  others. 

The  procedural  pattern  should  now  be  clear  for  synthesizing  an  nth 
order  system.   We  can  expect  to  have  n  initial  system  equations  in  2n 
unknowns.   There  will  be  (n-1)  necessary  conditions  to  satisfy,  and  (n-1) 
L-section  time  constants  to  be  determined  which  will  satisfy  these  condi- 
tions.  By  using  the  (n-1)  L-section  time  constants  and  redefining  some 
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of  the  cross  time  constants,  we  can  reduce  the  initial  system  equations 
to  n  equations  in  (n+1)  unknowns.   These  must  then  be  solved  for  n  un- 
knowns in  terms  of  the  remaining  one.   Assigning  a  value  to  this  remain- 
ing unknown  will  allow  us  to  go  back  and  solve  for  the  other  (2n-l)  un- 
knowns. 

It  should  be  noted  that  in  redefining  some  of  the  cross  time  con- 
stants to  derive  the  final  n  equations  in  (n+1)  unknowns,  it  is  neces- 
sary to  proceed  in  such  a  fashion  as  to  have  n  of  the  (n+1)  unknowns  be 
one  component  of  each  of  the  L-sections  and  the  other  unknown  to  be  the 
other  component  of  the  L-section  which  was  not  included  in  the  search 
for  allowable  time  constants.   As  an  example,  Eq.  (7.3a)  has  R„  and  C~ 
as  unknowns  since  we  found  allowable  values  for  the  first,  second  and 
fourth  L-section  time  constants,  and  Eqs.  (7.5)  contain  R  ,  R„,  R~ ,  and 

C,  as  unknowns. 
4 

Table  I  summarizes  the  observations  and  remarks  of  this  section. 
Table  I.   Characteristics  of  RC  Ladder  Networks 
of  the  First  Cauer  Form 


Network 

Basic 

#  of 

#  of 

#  of  Initial 

#  of  L-Section 

Order 

Criteria 

R's 

Cs 

System 

Time  Constants 

Conditions 

Equations 

to  be 

(Ineqs. ) 

Determined 

1 

0 

1 

1 

1 

0 

2 

1 

2 

2 

2 

1 

3 

2 

3 

3 

3 

2 

4 

3 

4 

4 

4 

3 

n 

n-1 

n 

n 

n 

n-1 

■» 
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9.    Computer -Aided  Network  Design. 

Use  of  a  digital  computer  greatly  facilitates  the  network  synthesis 
problem.   Once  a  program  has  been  written,  it  can  be  used  over  and  over 
again  to  obtain  any  number  of  suitable  networks  to  synthesize  a  given 
transfer  function.   It  can  provide  ranges  of  component  values  from 
which  appropriate  values  can  be  selected  to  synthesize  the  network. 
The  program  can  be  used  for  any  transfer  function  of  the  proper  form 
by  simply  changing  input  data  cards.   If  the  transfer  function  cannot 
be  physically  realized,  the  program  will  tell  that  too.   The  programs 
in  Appendices  II  and  III  illustrate  these  points. 

The  programs  in  Appendices  II  and  III  are  written  in  FORTRAN  63 
language  for  use  on  the  CDC  1604  digital  computer.   They  can  easily  be 
adapted  to  FORTRAN  60  and  FORTRAN  4  languages  with  just  minor  modifi- 
cations. 

The  program  in  Appendix  II  is  entitled  SYNTHNET.   It  is  written 
for  the  synthesis  of  a  third  order  RC  ladder  network  using  the  procedure 
of  section  6  of  this  paper.   Input  data  is  fed  to  the  computer  from  two 
data  cards  placed  at  the  end  of  the  card  deck.   Input  format  is  F10. 
Output  format  is  Ell.  4.   The  first  data  card  has  K  ,  K„ ,  and  K„  for 
the  given  transfer  function  punched  in  columns  1-10,  11-20,  and  21-30 
respectively.   Columns  31-40  contain  the  desired  initial  value  of  the 
L-section  time  constant,  a,    columns  41-50  the  amount  by  which  a   is  to  be 
incremented,  and  columns  51-60  the  desired  final  value  of  a.   The  second 
data  card  contains  the  desired  values  of  R  punched  in  each  set  of  ten 
consecutive  columns  from  column  1  to  column  80.   A  total  of  eight  values 
of  R  can  be  inputted  per  card.   Two  sets  of  component  values  will  be 
outputted  for  each  value  of  R  for  each  allowable  value  of  a  determined 
by  the  program  during  its  running.   The  corresponding  values  of  a   will 
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also  be  outputted.   Included  in  the  program  is  a  testing  procedure  to 
insure  that  the  component  values  which  are  calculated  actually  do 
synthesize  the  network.   The  original  K1 ,  K„,  and  K^  constants  are 
outputted  as  DK1 ,  DK2,  and  DK3 .   The  values  of  the  K' s  as  computed 
from  the  component  values  are  outputted  as  DCK1 ,  DCK2,  and  DCK3 . 
The  DK's  and  DCK's  should  be  the  same  respectively  if  the  component 
values  are  satisfactory. 

A  computer  run  was  made  with  the  following  as  input  data: 
KL   =   2.25 
K2  -   10.525 
K3   =   10.3 
Starting  point  of  a  =  0.05 
Step  size  of      a  =  0.05 
Stopping  point  of  a   =  10.0 
R   =   200  Kiloohms 
=   250  Kiloohms 
=     1  Megohm 
=   150  Kiloohms 
=  0.5  Ohms 
=   100  Megohms 
The  result  was  a  total  of  180  sets  of  component  values,  all  of  which 
will  synthesize  the  given  transfer  function. 

The  program  in  Appendix  III  is  entitled  DRAWCOMP.   It  is  written 
to  do  the  same  task  as  SYNTHNET,  except  that  its  output  is  in  the  form 
of  graphs.   Input  to  DRAWCOMP  is  the  same  as  for  SYNTHNET.   Each  com- 
puter run  will  provide  two  graphs,  one  for  resistors  and  one  for  capac- 
itors.  The  graphs  are  plotted  with  the  L-section  time  constant,  a, 
along  the  abscissa  and  the  component  values  along  the  ordinate.   We 
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can  determine  a  set  of  suitable  component  values  for  any  allowable 
value  of  a  by  extending  a  line  vertically  from  this  a  value  so  as  to 
intersect  all  of  the  curves.   The  proper  values  for  the  components  are 
then  found  at  the  intersections  of  the  vertical  line  and  the  various 
curves . 

A  computer  run  was  made  with  the  same  input  data  as  before,  ex- 
cept that  the  step  size  of  a  was  changed  to  .0004  and  the  starting  and 
stopping  points  changed  to  0.85  and  1.55  respectively,  since  these  are 
limiting  values  of  a  as  shown  by  the  output  data  of  SYNTHNET.   Results 
of  the  run  are  shown  in  Figs.  5. a.  through  10. b.   It  is  noted  that  all 
of  the  curves  have  the  same  shape,  as  might  be  expected.   In  fact,  if 
we  superimpose  the  curves  of  Figs.  9. a.  and  10. a  for  the  capacitors 
and  Figs.  9.b  and  10. b  for  the  resistors  it  will  be  seen  that  they  are 
identical.   The  only  thing  different  is  the  ordinate  scale  factor.   Van 
Valkenburg   states  in  his  text  that  if  each  resistor  in  the  network  is 
multiplied  by  some  constant,  d,  and  each  capacitor  by  the  reciprocal 
of  that  same  constant,  1/d,  then  the  driving  point  functions  of  the  net- 
work will  have  been  scaled  by  that  factor  but  the  voltage-ratio  and 
current-ratio  transfer  functions  are  not  affected.   This  fact  is  borne 
out  in  the  graphs.   If  we  input  one  desired  value  of  R1  to  the  computer, 
we  will  get  a  set  of  curves.   If  we  now  input  another  value  of  R1 ,  then 
in  effect  we  will  have  scaled  the  original  R1 .   The  second  set  of  curves 
will  give  the  same  component  values  as  the  first  set  of  curves,  but  they 
will  be  scaled.   It  can  be  stated  then  that  a  single  set  of  curves  for 
a  given  transfer  function  will  provide  a  complete  set  of  data  for 


M.  E.  Van Valkenburg,  Introduction  to  Modern  Network  Synthesis 
(New  York:   J.  Wiley  and  Sons,  1960),  pp.  48-53. 
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synthesizing  that  transfer  function.   Inputting  additional  values  for 
R   to  the  program  merely  provides  one  with  a  convenient  means  of 
scaling  the  entire  set  of  network  components. 
10.   Conclusions . 

Up  to  this  point  we  have  considered  ladder  networks  made  up 
entirely  of  RC  elements.   The  principles  of  duality  can  be  applied  if 
we  wish  to  consider  RL  ladder  networks  too.   The  procedure  outlined  can 
be  used  equally  well  for   synthesizing  RL  networks  from  a  given  current- 
ratio  transfer  function  by  replacing  all  resistances  of  the  RC  network 
by  inductors  and  all  capacitors  by  conductances. 

It  has  been  assumed  thus  far  that  any  impedance  connected  to  the 
output  terminals  of  the  network  will  be  of  such  magnitude  as  to  have 
negligible  loading  effect  on  the  network.   The  actual  effect  of  various 
types  of  load  impedances  upon  the  currents  and  voltages  within  the  net- 
work can  be  readily  determined  by  using  the  results  of  section  4  and 
continuants. 

It  has  been  shown  that  there  are  a  large  number  of  sets  of  compo- 
nent values  which  can  be  used  to  synthesize  a  ladder  network  for  a 
given  transfer  function.   Each  of  these  sets  will  have  a  different  in- 
put impedance.   Eq.  (4.2)  could  be  included  in  a  computer  program  for 
network  synthesis  to  give  the  input  impedance  for  each  set  of  component 
values  determined.   The  scaling  scheme  previously  discussed  could  also 
be  included  in  the  program  as  desired. 

A  logical  extension  of  the  work  outlined  in  this  paper  would  seem 
to  be  the  writing  of  computer  programs  to  synthesize  higher  order  sys- 
tems, investigate  the  loading  effect  on  the  networks  for  various  net- 
work terminations,  determine  input  and  output  impedances  of  the  net- 
works for  each  set  of  component  values,  and  perhaps  develop  and  program 
the  equations  for  synthesizing  RL  networks. 
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APPENDIX  I 


THEOREMS  OF  CONTINUANTS 


The  following  theorems  are  taken  from  the  text  by  Bartlett 

l"   K^a15a_,,.a  )   =  K  ( a  ,a   , , . . .a„ ,a  j 
1   2    n         n  n~l     2   1 

2.   k(a, ,a_, . „ ,a  )  =  a.k(a0, . . .a  )  +  k(a„,...a  ) 
12'     n      1    2     n         3     n 

3„   k(a1}...a  )   =  k(a.,...a  )k(a  .....a  ) 
1     n         1     P    P+l    n 

+  k(a.,...a   ,)k(a   ...a  ) 
1     p-1    p+2    n 


4. 


k(a.,...a  )k(a~,...a   .. )   =  k(a.,...a   n)k(aOJ...a  ) 
1     n    2     n-1         1     n-1    2     n 


+  (-D 


k(a1je..a  )   =  k(a1s.„.a   ,)k(a  ,„..a  ) 
v  1'    n         1'    r-1    r'    n 


+k(a15...a   0)k(a  ,,,... a  ) 
1     r-2    r+1     n 


6. 


k(0,a_,...a  )   =  k(aOJ...a  ) 
2     n         3     n 


A.  C.  Bartlett. ,  The  Theory  of  Electrical  Artificial  Lines 
and  Filters  (New  York:   J.  Wiley  and  Sons,  1931),  p.  47 
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APPENDIX  U 


2  10  67  PAGE  NO. 


PROGRAM  SYNTHNET 

DIMENSION  D»('0 ) ,R2(2) ,C2(2) 

DO  999  Jsl,50  

999  DR( J>=0 . 0 

DO  998  K=l,2 

R2(K)=0.0 
998  C2(K):0. 0 

READ  5'0,DlO.r)K2.D<3.STARTA,STEPA.ST0PA 
50  roRMAT(6rl0)  .  _ 

READ  60,DR(1),CR(2),DR(3),DR(*)»DR(5),DR(6)»DR(7),DR(8) 

60  roRMATono) 

A*STARTA 

100  Y*DK2/a-a-(2. #CK1 )/( A**2i 
IFlY)  101,1 01,102 

101  A=A*STEPA 
IFISTOPA-A*  125,125,100 

102  X«DK3-?.«a.Dk1/ ( a**2) 
IF  CX  )  101,101,103 

1 03  Tl3=x-y 
IF(Tl3)101,lni.lO< 

10*  W*( Y**2)-(4..Dk1«t13)/( A**2) 
IF(W)  101,105,105 

105  CONTINUE 

DO  H6  Isl,8 
IF(DR( I ) )125, 101,115 
115  Rl=DR( I ) 
C3=T13/Rl 
Ct=A/Rl 
R3=A/C3 

r2(1);(y*sq^tF'^))/(2,.c3) 
C2( l)rDKl/' (a»*2)*r2(  1)  ) 
R2(2)=(Y-SQRTF(W))/'2.*C3) 
IF(R2(2)  )107. 109,106 

106  c2(2)rDKl/( ( A*«2)«p2(2) ) 
GO  TO  111 

107  PRINT  me 

108  FORMAT l  ///,  10<, 17HR2( 2)  IS  NEGATIVE) 
GO  TO  113 

109  PRINT  110 

110  FORMAT  t ///,  10x ,13hb2(2)  js  7ERO> 
GO  TO  113 

111  PRINT  112,a,k1,r2(2) ,R3,cl,c2f 2) ,c3 

112  FORMAT  <///,10m,2hai,E11.4,1P*,3hb1=,E11.4.10X,6Hr2(2)::,e11.4, 
H0X,3HK3r,pll.4,/,33x,3HCl=.f:ll.4,10x,6HC2(2)r,Ell.4,10x,3Hc3r, 
2F11.4) 

DCK1iR1*C1*-j2(2)«c2(2).r3»c3 

nCK2  =  Fl*cl*'-'2(  2  ).Q2r  2  )*Rl.cl*o2<  2)  *c3*r1*c1*r3#c3  +  r1.c2(  2)  «r3*c3 
l+«2( 2 i .q2( ? ) ,^3*c3 
ncK3;ul.cl*Rl»c2(2)*Rl*C3*R2(2)«c2(2)*R2(2)*c3+R3,c3 
PRINT  120,  nCKl . DC<2 . DCK3 

120  F0WMAT(/>33x,5h0C<l-.Ell.4.8x.5HDCK2=,Ell.4,llx,5HDC'<3r,Ell.4) 
PRINT  121 , n«l . CK2, dk3 

121  FORMaT(33x,4mo«<1=.E11.4,9X,4HDK2:,e11.4>12x,4hdk3=,e11.<») 

113  print  114, a,r1.r2(1]  ,r3,c1,c?( 1) ,C3 

114  FOHMAT  (/,inx,2HA»,E1l.4,lfV.3MRlz,Ell.4,10x,6HR2(l)r,Ell.4. 
110x.3hk3=,p11.4,/,33x,3hc1=.f11.4,10x,6hc2(1):,fH.4,10x,3hc3=, 
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1  1* 


1  =  P1*C1*^?(1  i  * C 2 1  1  )  * «} 3 «  c 3 

?=frl*Cl«<''(l)«C'?,1i*^1«C1*t?2(1)*C:,*^l*C1*^3*c3  +  Rl*c2(l)*p3*c3 
(1  ) . c  ? ( 5  i .  j  '  »  r  5 
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DCK?:     1.0525F     01 
PK2=     1. 0525E    01 

p2( 1)=  7 . T425t  06 
C2(  1  ) z     4 . 3601F-07 

DC*2=  1.0525p  01 
DK2r  1.0525F  01 


p3;  1.1686F  06 
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nCK^s  1  .  030  0b  01 
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p3=  1.1686F  06 
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DCK2z  1.0525e  01 
DK2:  1.0525E  01 


r3s  1.6574?  05 
c3r  6.3350E-06 

OCK3s  1.03006  01 
DK3:  1.0300e  01 

qi=    1.6574f  05 
n3=  6.33506-06 

DCK3s  1.0300E  01 
DK3s  1.Q3006  01 


as  1. 05Q0F  on 


Bis  2.5po06  05 
Cl:  4.20006-06 

DCKls  2.?5no=  00 
DKls  2.?5006  00 


r2(2):  1.18956  05 
C?<?>=  1.71576-05 

DCK2s  1.0525e  01 
DK2r  1.0525E  01 


r3s  2.0718F  05 
C3s  5. 06806-06 

DCK3s  1.0300F  01 
DK3:  1.q3oOp  01 
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At  1.0500E  00  Rls  2.5000E  05  R2<1).  8.«636E  05  r3.  2.0718F  05 

Cl=  <.2000E-06  C2<1)«  2,4113e-06  C3»  3.0680E-06 


DCKls  2.2500E  "00         DCk2»  1.0525e  01  DCK3s  1.0300E  01 

DKlr  2.2500E  00  DK2s  1.05?5e  01  DK3  =  1.Q300E  01 


A=  1.0500E  00  Rls  1.0JQQE  H6     R2(2)8  3^7579B    05 -g3g  ft.2g72E  05 

Cls  1.0500E-06  C2<2>=  «.289*E-06  C3»  1.2670E-06 

DCKls  2.25QOE  00  DCK2s  l.Q525E  01 BCH3^  i.  fl3oOE  01 

DKlr  2.2500E  00  DK2=  1.0525E  01  0*3  =  1.0300E  01 

A=  1.05QOE  00           Rls  1.0000E  Q6- R2(D«  3.3854e  06  R_3«J.2a72E  05 

Clt  1.0500E-06  C2<D=  6.02g2E-07  C3«  1.2670E-06 

DCKl:  2.2500E  QQ  _.  _  JJCK_2s_l.  052$£  01  DC«3=  l.fl3oOE  01 

DKlr  2.2500E  00  DK2=  1.0525E  01  DK3r  1.Q300E  01 


Ar  1.0500E  00  Rlr  1.5000E  05  R2(2)r  7.1368E  0*  r3«  1.243lE  05 

Clr  7.0Q0QE-Q6  j;2{2)  =  2.8596E-05  C3*  8..««67E-06 

DCKlr  2.2500E  00  DCK2s  1.0525E  01  DCK3=  1.0300F  01 

DKl=  2.2500E  00  DK2r  1.0525E  01  DK3r  1.03OOE  01 

Ar  1.0500E  00  Plr  1.5000E  05  R2(l)r  5.0781E  05  r3»  1.2431e  05 

Cl=  7.0000E-06  C2(l)r  4,oiB8E-06  C3»  8.4467E-06 

DCKlr  2.2500E  00  DCK2r  1.0525E  01  DCK3r  1.0300E  01 

DKl=  2.25Q0E  00  DK2r  1.0525E  01  DK3=  1.0300E  01 

Ar  1.0500E  00  Plr  5.0000E-01  r2<2)=  2.3789F-01  r3«  4.1436E-01 

Clr  2.1000E  00  C2(2)r  8.5787E  00  C3*  2.5340E  00 

DCKlr  2.2500E  00  DCK2r  1.0525E  01  DCK3r  1.0300E  01 

DKlr  2.2500E  00  DK2=  1.0525E  01  DK3r  1.0300E  01 

Ar  1.0500E  00  Rlr  5.0000E-01  R2(l)r  1.6927E  00  r3«  4.1436E-01 

Clr  2.1000E  00  C2<l)r  1.2Q56E  00  C3»  2.5340E  00 

DCKlr  2.2500E  00  DCK2r  1.0525E  01  DCK3r  1.0300E  01 

DKlr  2.2500E  00  DK2=  1.0525E  01  DK3r  1.Q300E  01 


Ar  1.0500E  00  Rlr  1.0000E  08  R2(2)=  4.7579E  07  r3«  8.2872e  07 

Cl=  1.0500E-08  C2(2)S  4.2894E-0*  C3*  l.?670F-08 

DCKlr.  2.2500E  00  DCK2s  1.0525E  01  DCK3r  1.0300E  01 

DKlr  2.2500E  00  DK2r  1.0525E  01  DK3r  1.0300E  01 

Ar  1.0500E  00  Rls  1.0000E  08  R2(l)r  3,3854E  08  r3«  8.2fl72E  07 

Clr  1.0500E-08  C2(l)r  6.0282E-09  C3»  1.2670F-08 

DCKlr  2.25Q0E  00  DCK2r  1.0525e  01  DCK3r  1.0300E  01 

DKlr  2.2500E  00  DK2r  1.0525E  01  DK3r  1.q300e  01 

Ar  1.1000F  00  Rlr  2.0000E  05  R2(2)s  9.1435E  04  p3s  1.4752e  05 

rlr  5.5nnnp-06  r"2(2)s  2.n337p-o5  r3s  7.4566F-06 
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As  1.10Q0E  00 


DCKls  2.2500E  00 
DKls  2.25001?  00 

Rls  2.0000E  05 
Cls  5.5000E-06 

DCKls  2.25006  00 
DKls  2. 25006  00 


DCK2s  1.0525E  01 
DK2=  1.05256  01 

R2<l)s  5.4547E  05 
C2(l)=  3.4090E-06 

DCK2s  1.0525E  01 
DK2r  1.05?5e  01 


DCK3*  1.0300E  01 
dk3=  1.0300E  01 

r3«  1.4752f  05 
c3»  7.45666-06 

DCK^r  1.0300E  01 
DK3r  1.03Q0E  01 


A=  1.1000E  00 


A=  1.1000F  00 


5l:  2.5000E  05 
Cl=  4.4000F-06 

DCKls  2.25006  00 
DKls  2.2500E  00 

Bis  2.5000E  05 
Cls  4.4O00E-06 

DCKls  2.25006  00 
pKls  2.25006  00 


R2(2)S  1.1429E  05 
C2(2)s  1.6270E-05 

DCK2=  1.0525e  01 
DK2=  1.o5?5e  01 

R2(l)=  6.8184E  05 
C2(l)=  2.72725-06 

DCK2=  1.0525E  01 
DK2s  1.05256  01 


r3i  1.844of  05 
(;3»  5.9653f-06 

DCK3=  1.0300E  01 
DK3i  1.0300F  01 

r3»  1.844of  05 
C3»  5.9653f-06 

DCK3r  1.0300F  01 
DK3z  1.Q300F  01 


A:  1.10Q0F  00 


a=  l. iooof  on 


Rls  1.00006  06 
Cls  1.10006-06 

DCKl=  2.25005  00 
DKlr  2.25006  00 

Bis  1.00006  06 
Cls  1.1000F-06 

DCKls  2.25006  00 
flKls  2.25006  00 


r2(2)s  4.57186  05 
C2(2)=  4.0674E-06 

DCK2=  1.05256  01 
DK2s  1.0525E  01 

r2(1)=  2.7?746  06 
C2(l)=  6.81806-07 

DCK2=  1.05256  01 
DK2=  1.05256  01 


r3=  7.37606  05 
C3=  1.4913f-06 

DCK3s  1.03006  01 
DK3s  1.0300e  01 

r3s  7.3760F  05 
c3*  1.4qi3f-06 

DCK3s  1.03006  01 
DK3r  1.03Q06  01 


as  1.  10006  00 


as  1.  1000F  00 


pis  1.50006  05 
Cls  7.3333e-06 

[)CKls  2.2500E  00 
DKli  2.25006  00 

Bis  1.5000E  05 
Cl=  7.3333e-06 

DCKls  2.2500E  00 
DKls  2.?500t  00 


r2(2)=  6.8576E  04 
C2(2)s  2. 7U66-05 

DCK2r  1.05256  01 
DK2s  1.05?5e  01 

*2<  1)=  4.09106  05 
C2(l)=  4.S453E-06 

DCK2s  1.0525f  01 
DK2r  1.05256  01 


r3s  1.1064f  05 
r;3=  9.942lf.06 

DCK3r  1.03006  01 
r)K3=  1.03006  01 

r3s  1.1064F  05 
c3=  9.9421F-06 

DCK3s  1.0300F  01 
Dk3=  1.03q06  01 


As  1. 10006  00 


as  l.inonp  0" 


Rls  5.  OOOOF-01 
Cls  2.20006  on 

DCKls  2.25006  00 
DKls  2.2500F  00 

Rls  5.onnoF-oi 

rls  2.2000P  00 

DCKls  2.25005  00 
DKls  2.P5006  00 


r2(2)s  2.28596-01 
C?( 2)s  8.  134qf  00 

DCK2s  1.05255  ol 
DK2r  1.0525F  01 

R2(l)s  1.3fc37F  00 
C2(l)r  1.36366  00 

Dck2=  1.0525g  01 
DK2s  1.05256  01 


r3s  3.6880F-01 
CJs  2.9g26F  00 

DCK3=  1.03006  01 
DK3s  1.03006  01 

r3=  3.68806-01 
f^3s  2.9826F  00 

DCK3s  1.03006  01 
nK3z  1.03006  01 
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A:  1.1000E  DO 


A=  l.lOOOE  00 


Rls  1.0QD0E  dfl  _ 
Cis  1.1O00E-08- 

DCKls  2.250QE_DQ_ 
DKls  2.25nOE  00 

Rls  1.0000E  08 
Cls  1.1000E-08 

DCKls  2.2500E  00 
DKls  2.2500E  00 


j2(2)=  <;57lBE  07_ 
C2(2)s  ".0674E-08 

DCK2»  1.Q525E  01 
DK2r  1.0525E  01 

R2(l)r  2.7274fc  08 
C2(l)=  6.8180E-09 

DCK2=  1.0525e  01 
DK2r  1.0525E  01 


J?3«  7.3760E  07 
C3r  1.49i3e-08 

DC*3-*  1.0300E  01 
OK3i  1.0300E  01 

r3j  7.3760E  07 
c3i  1.49l3p-08 

BCK3i  1.0300E  01 
qk3=  1.0300F  01 


A:  1.1500E  00 


As  1.15Q0E  00 


Rl=  2.0000E  0' 
Cls  5. 7500E--OA  

DCKls  2.2500E  00 
DKls  2.2500E  00 

Rl=  2.0000E  05 
Cls  5.7500E-06 

DCKls  2.25Q0E  00 
DKls  2.2500E  00 


R2(2)3  8.B*18E  04 
C2(2x*  1.92<2e-05 

DCK2x  1.0525E  01 
DK2:  1.0525E  01 

R2(l)=  4.5297E  05 
C?(l)s  3.7559E-06 

DCK2s  1.0525E  01 
DK2r  1.Q525E  01 


r3«  1.3536F  05 
Cil  8.4957E-06 

DCK3s  1.0300F  01 
DK3s  1.Q3Q0E  01 

r3»  1.3536F  05 
C3s  8.^9576-06 

DCK3=  1.0300E  01 
DK3=  1.0300E  01 


as  1.1500E  00 


As  1.1500E  00 


His  2.5000E  05 
Cls  4.6000E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Rls  2.5000E  05 

Cls  4.6000F-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 


R2(2)s  1.1052E  05 
C2(2)=  1.53p3E-05 

DCK2s  1.0525F  01 
DK2s  1.05256  01 

R2(l)=  5.6622E  05 
C2(l)s  3.  00*7E-06 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 


r3s  1.6920F  05 
c3»  6.7966E-O6 

DCK3s  1.0300E  01 
DK3i  1.0300E  01 

r3»  1.6920F  05 
C3s  6.7966E-06 

OCK3s  1.0300E  01 
DK3:  1.0300E  01 


As  1.15Q0 E  00 


As  1.1500E  00 


Rls  1.0000F  06 
Cls  1.1500E-06 

DCKls  2.2500E  00 
DKls  2.25n0E  00 

Rl=  1.0000E  06 
Cls  1.1500E-06 

DCKls  2.2500;  00 
DKls  2.2500E  00 


R2(2)s  4.4209E  05 
C2(2)s  3.8483f-0^ 

DCK2s  1.0525E  01 
DK2s  1.05256  01 

R2(l)s  2.2649e  06 
C2(l)s  7.5118E-07 

DCK2s  1.0525P  01 
DK2s  1.0525E  01 


r3«  6.7681E  05 
c3«  1.699le-06 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 

r3=  6.7681E  05 
C3s  1.6991F-06 

DCK3s  1.0300E  01 
DK3s  1.Q300E  01 


As  1.1500E  00 


as  1.  15Q0E  00 


ols  1.5000E  05 
Cls  7.6667E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

els  1 . 5000E  05 
rl:  7.6(SA7P-06 


R2(2)r  6.6314E  04 
C2(2)=  2.5656E-05 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 

r2(1>=  3.3973e  05 
r2(  1  )s  5  .  (107OF-06 


r3=  1.0152f  05 
C3S  1.1328E-05 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 

r3.  1.0152e  05 
r3s  1.132AF-05 
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ZX 


DCKls  2.2500E  00 
DKls  2.2500E  00 


DCK2S  1.0525e  01 
DK2:  1.05P5E  01 


ncK3s  1.0300E  01 
DK3:  1.0300E  01 


A=  1.1500E  00 


A  =  1.15Q0F  00 


pis  5. OOOOE-Ol 

c-l-s  2.3nonp  oo 

DCKls  2.2500E  00 
DKls  2.?500E  00 

Rls  5. 0000E-01 
Cl=  2.3000R  00 

DCKls  2.25Q0E  00 
DKl=  2.2500E  00 


R?(2)S  2.2105E-01 

C2(2)s  7.6967E  00 

DCK2«  1.0525E  01 
DK2s  1.0525E  01 

R2(l)s  1.1324E  00 
C2(l)=  1.5o24E  00 

DCK2s  1.0525g  01 
DK2=  1.0525E  01 


r3«  3.384of-01 
c3a  3.3983f  00 

DCK3s  1.0300E  01 
DK3s  1.0300F  01 

r3.  3.384oe-01 
c3s  3.39g3F  00 

DCK3=  l,0300E  01 
DK3s  1.0300E  01 


as  1.15Q0F  00 


As  1.15Q0E  00 


Kir  1.  0O0OE  08 
Cls  1.1500E-08 

DCKls  2.2500E  00 
DKls  2.2500E  00 

His  1.0000E  08 
Cls  1.1S00E-08 

DCKls  2.2500E  00 
DKls  2.2500E  00 


r2(2)s  4.4209E  07 
C2(2)=  3.b483E-OB 

DCK2s  1.0525E  ol 
DK2s  1.0525E  01 

R2(l)s  2.26«9E  08 
C2( 1 ) s  7. 5118F-09 

DCK2s  1.0525E  01 
DK2r  1.0525E  01 


r3=  6.7681E  07 
C3s  I.699IF-OB 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 

r3b  6.7681E  07 
c3=  I.699IF-O8 

DCK3s  1.0300E  01 
DK3s  1.03Q0E  01 


as  1.2000F  00 


Ar  1.2000E  00 


uls  2.0000E  05 
Cls  6.0000E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

pis  2.0000F  05 
Cl=  6.0000F-06 

DCKls  2.2500E  00 
DKls  2.2500F  00 


R2(2)=  8.60<0fc  04 
C2(2)S  1.8160F-05 

DCK2s  1.0525E  01 
DK2=  1.0525E  01 

R2( l)s  3.8400F  05 
C2(l)s  4.0690E-06 

DCK2=  1.0525F  01 
DK2s  1.0525F  01 


r3=  l.?687F  05 
c3:  9i45e3E-06 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 

r3=  1.2687F  05 
c33  9.4583f-06 

DCK^s  1.0300E  01 
DK3s  1.Q300F  01 


as  1.2000F  00 


«s  1.2000F  00 


-ilz    2  .  5000F  05 
Cls  4,9000E-06 

DCKls  2.2500E  00 
DKls  ?.?500F  00 

ol:  2. 5000F  05 
Cls  4.8000E-06 

DCKls  2.2500E  00 
DKls  2.?500E  00 


R2(2)s  1.0755E  05 
C2(2)=  1.4528E-05 

DCK2s  1.05?5e  01 
DK2=  1.0525E  01 

R2(l)s  4.8001E  05 
C2(l)=  3.2552E-06 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 


r3s  1.5a59F  05 
C3r  7.5667F-06 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 

r35  1.5859F  05 
t-3s  7.5667F-06 

DCK3s  1.0300F  01 
DK3s  1.03Q0E  01 


as  1.2000E  00 


pis  1.0000E  0* 
Cls  1  .  2n0OE-O6 

nr*rlr    2    25nn:    nn 


r2(2)=    4.302oe    05 
C2(2)=    3.A320E-06 

nr-K2s    1  .  n5?5n    nl 


r3s  6.3436P  05 
c3=  1.8917F-06 

n  r  if  3 1  s  1 .  n  3  n  n  p  01 


is- 


DKl=  2.2500E  00 


DK2=  1.0525E  01 


DK3:  1.0300E  01 


A=  1.2000E  00 


Rl:  1.0000E  06 
Cls  1.2000F-06 

DCKls  2.250QE  00 
DKls  2.2500E  00 


_fi2(l)r  l.g200E  06 
C2(l)=  8.1379E-07 


_OCK2=  1.0525e  01 
DK2:  1.o5?5e  01 


r3s  6.3436E  05 
C3S  1.8917E-06 

DCK3r  1.0300E  01 
DK'=  1.0300E  01 


A  =  1.2000E  00 


A  =  1.2000E  00 


Pis  1.5000E  05 
Cls  8.0000E-06 

OCKls  2.25Q0E  00 
DKls  2.2500E  00 

Rls  1 .  50001=  05 

Cls  8.0000E-Q6 

DCKls  2.2500E  00 
DKls  2.P5Q0E  00 


R2(2)S  6.4530e  o* 
C2(2)r  2.4214E.05 

DCK2s  1.0525e  01 
0*2:  1.0525E  01 

R2(l)s  2.8800E  05 
£2(l)s  5.4253E-06 

DCK2s  1.0525E  01 
DK2:  1.0525E  01 


r33  9.5l54E  o4 
C3s  1.2611F-05 

DCK3s  1.0300E  01 
OK3=  1.0300E  01 

r3,  9.5154E  0* 
C3«  1.2611F-05 

DCK3r  1.0300E  01 
DK3r  1.Q300E  01 


As  1.2000E  00 


as  1.2000E  00 


Rls  5.0000E-01 
Cls  2.4000E  00 

DCKls  2.2500E  00 
DKls  2.2500F.  00 

ols  5  .  000flP-0l 
Cls  2.4000E  00 

DCKls  2.2500E  00 
DKls  2.2500E  00 


r2(2)=  2.1510E-01 
C2(2)S  7.2641E  00 

DCK2s  1.0525e  01 
DK2=  1.0525E  01 

R2(l)=  9.6001E-01 
C2(l)s  1.6276E  00 

DCK2:  1.0525E  01 
DK2:  1.0525E  01 


r3=  3.1718E-01 
c32  3.7833E  00 

DCK3s  1.0300E  01 
DK3r  1.0300E  01 

r3=  3.1718E-01 
c3=  3.7833e  00 

DCK3r  1.0300E  01 
nK3=  1.0300F  01 


As  1.2000E  00 


As  1.2000F  00 


Rls  1.0000E  08 
Cls  1.2otlOfc-08 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Rls  1.0000E  08 
Cls  1.2000P-08 

DCKls  2.2500=  00 
DKls  2.2500E  00 


R2(2)=  4.302oe  07 
C2(2)=  3.632oe-08 

DCK2r  1.0525E  01 
DK2s  1.0525E  01 

p2(l)=  1.9200E  08 
C2(l)=  8.1379E-09 

DCK2:  1.0525c  ol 
Dk2=  1.0525E  01 


r3s  6.3436E  07 
C3s  1.8917E-08 

DCK3s  1.0300E  01 
DK3r  1.Q300E  01 

R3a  6.3436E  07 
C3s  1.8917F-08 

DCK3s  1.0300F  01 
DK3r  1.0300E  01 


As  1.25Q0F  00 


as  1.2500F  00 


Pis  2  .  0000P  05 
Cls  6.2500E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Rls  2.0000E  05 
Cls  6.2500E-06 

DCKls  2.2500F  00 
DKls  2.25Q0E  00 


R2(2)=  8.4263E  04 
C2(2)z  1.7089E-05 

DCK2r  1.0525E  01 
DK2=  1.0525E  01 

R2(l)r  3.3o23e  05 
C2(l)s  4.3606E-06 

DCK2=  1.0525F  01 
DK2r  1.0525E  01 


r3s  1.2077E  05 
C3s  1.0350F-05 

DCK3:  1.0300E  01 
DK3s  1.0300E  01 

R3r  1.2077E  05 
c3s  1.0350F-05 

DCK^s  1  .  0300E  01 
DK3=  1.Q3Q0F  01 
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As  1.25Q0E  00 


A=  1.2500E  00 


Rls  2.5000E  05 
Cl=  5.0000E-06 

DCKls  2.25006  00 
DKl=  2.2500E  00 

pi*  2.5000E  05 
Cl:  5.0000E-06 

DCKl=  2.2500E  00 
DKls  2.2500E  00 


R2(2)S  1.0533e  05 
C2(2)s  1.3672E-05 

DCK2»  1.0525E  01 
DK2r  1.0525E  01 

R2(l)=  4.1279E  05 
C2(1)S  3.4885e-06 

DCK2=  1.0525E  01 
DK2:  1.0525E  01 


r3i  1.5097E  05 
c3x  8.2800E-06 

DCk3=  1.0300E  01 
OK'r  1.0300E  01 

p3«  1.5097E  05 
c3»  8.2flOOE-06 

DCk3=  1.0300F  01 
DK3=  1.0300E  01 


As  1.25Q0E  00 


as  1.25Q0E  00 


Sis  l.OOOOE  06 
Cls  1.2500E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Pis  l.OOOOE  06 
Cls  1.2500E-06 

PCKls  2.2500E  00 
DKls  2.P500E  00 


r2(2)=  4,2131E  05 
C2(2)=  3.4179E-06 

DCK2s  1.0525E  01 
DK2r  1.0525E  01 

R2(l)s  1.A512E  06 
C2(l)s  8.72l2e-07 

DCK^s  1.0525g  01 
DK2s  1.0525F  01 


r3s  6.0386E  05 
c3»  2.0700F-06 

DCK3s  1.0300E  01 
OK3s  1.0300E  01 

r3s  6.0386E  05 
c3a  2.0700F-06 

DCK3s  1.0300E  01 
nx3s  1.03Q0E  01 


as  1.2bonE  00 


as  1.2500E  00 


Rls  1.5000E  05 
Cls  8.3333e-06 

DCKls  2.2500=  00 
DKls  2.?500E  00 

Rls  1.5n00E  05 
Cls  8.3333E-06 

DCKls  2.2500=  00 
DKls  2.2500E  00 


R2(2)=  6.3197E  04 
C2(2)=  2.2786E-05 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 

R2(l)r  2.4767E  05 
C2(l)s  5.8141E-06 

DCK2s  1.05?5p  01 
DK2s  1.0525E  01 


r3s  9.0580E  04 
p3=  1.3fl00E-05 

0ck3s  1.0300E  01 
DK3=  1.0300E  01 

r3i  9.0580F  04 
C3=  1.3800F-05 

DCK3s  1.0300F  01 
DK3s  1.Q3Q0E  01 


As  1.2500F  00 


a=  1.2500E  00 


Ml=  5.0000F-01 
Cls  2.5000F  00 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Jl:  5.  OOO0F-01 
Cls  2.5000E  00 

DCKls  2.2500E  00 
DKls  2.2500E  00 


R2(2)s  2.1066E-01 
C2(2)r  6.fl358E  Of) 

DCK2s  1.0525E  01 
DK2r  1.0525E  01 

r2(1)s  8.2558E-01 
C2(l)s  1.7442F  00 

DCK?=  1.0525E  01 
DK2r  1.0525E  01 


r3-  3.0193f-01 
c3r  4.1400E  00 

DCK3s  1.0300E  01 
DK3z  1.0300E  01 

r3s  3.0l93p-01 
c3s  4.1400E  00 

HCK3s  1.0300F  01 
DK3r  1.03Q0E  01 


as  1.2500E  00 


as  1.25Q0E  00 


Rls  1. 0O00E  08 
Cls  1.2500F-08 

DCKls  2.2500E  00 
DKls  2,?5n0F  00 

pis  1.  0000c.  08 
cls  1.2500E-08 

nr<l:  2  25nnc  nn 


R2(2)s  4.2131E  07 
C2(2)=  3.4i79p-08 

DCK2s  1.0525e  01 
DK2r  1.0525E  01 

R?(l)=  1.6512E  08 
C2< l)s  g. 72l2E-0p 

nrK2s  l.n525c  nl 


r3s  6.0386F  07 
C3»  2.0700F-08 

DCK3s  1.0300E  01 
D<3-  1.0300E  01 

r3=  6.03e6p  07 
C3s  2.0700F-08 

nrK3:  l .  n 3 n  n  =  n i 
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DKls  2.2500E  00 


DK2z  1.0525e  01 


0K3s  l.O^OOp  01 


A:  1.3000E  00 


A=  1.3000F  00 


Rli  2.0000E  05 
Clj  6.5000E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Rls  2.0000E  05 
Cl=  6.5000E-06 

DCKlr  2.2500E  00 
DKls  2.2500E  00 


R2(2)s  8.3oe3E  04 
_C2(2)=  l^fl?4E-o5 

DC»<2s  1.0525e  ol 
DK2s  1.0523E  01 

R2(l)=  2.B676E  05 
C2(l)=  4.6427E-06 

DCK2r  1.0525E  01 
DK2s  1.0525E  01 


r3,  1.1632f  05 
c3=  1.1176E-05 

DCK3s  1.0300E  01 
0«3r  1.03Q0F  01 

r3s  1.1632e  05 
C_3,  1.H76E-05 

DCK3=  1.0300E  01 
0X3=  1.Q300E  01 


A:  1.3000E  00 


A=  1.3000E  00 


Rl=  2.5000E  05 
Cls  5.2000E-06 

DCKls  2.25Q0E  00 
DKls  2.25Q0E  00 

Rls  2.5000E  05 
Cls  5.2000E-06 

DCKls  2.2500E  00 
DKls  2.25Q0E  00 


R2(2)s  l.n3§5E  05 
C2(2)=  1.2820E-05 

DCK2s  1.0525g  ol 
0K2s  1.0525E  01 

~r2(1)=  3.58<6e  05 
C2(l)=  3.71<2e-06 

DCK2=  1.o525e  01 
DK2s  1.05?5e  01 


r3=  1.4540E  05 
C3«  8.9«08F-06 

DCK3r  1.0300E  01 
qK's  t.0300E  01 

r3=  1.454oe  05 
C3«  8.9408E-06 

nCK3=  1.0300E  01 
OK3r  1.Q300E  01 


A=  1.3000E  00 


A=  1.3000E  00 


Rls  1.0000E  06 
Cls  1.3000E-06 

DCKls  2.2500E  00 
DKls  2.?500E  00 

Rls  1.0000E  06 
Cls  1.3000E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 


r2(2,,  4.1542E  05 
C2(2)s  3.2o*9E-06 

DCK2s  1.0525e  01 
DK2s  1.0525E  01 

r2(1)-  1.4338E  06 
C2(l)s  9.?854E-07 

DCK2s  1.0525e  01 
0K2s  1.0525E  01 


r3,  5.8160E  05 
C3S  2.2352p-06 

DCK3s  1.0300E  01 
nK3s  1.0300E  01 

r3i  5.8160E  05 
c3«  2.2352e-06 

DCK3s  1.0300E  01 
DK3s  1.03QOE  01 


A  =  1.3000F  00 


As  1.3000E  OP 


Pis  1.5000E  05 
Cls  8.6667E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Rls  1.5000E  05 
Cls  8.6667E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 


r2<2)s  6.2312E  04 
C2(2)=  2.1366E-05 

DCK2s  1.0525e  01 
DK2s  1.o5?5e  01 

R2(l)s  2.1507E  05 
C2(l)s  6.1903E-06 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 


r3s  8.724oe  04 
C3«  1.49oiE-05 

DCK3s  1.0300E  01 
DK3s  1.Q300E  01 

r33  8.724of  04 
C3«  1.4901F-05 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 


As  1.3000E  00 


Rls  5.0000E-01 
Cls  2.6000E  00 

DCKls  2.2500E  00 
DKls  2.2500E  00 


r2(2)=  2.0771E-01 
C2(2)s  6.4098E  00 

DCK2s  1.0525E  01 
DK2=  1.0525E  01 


r3»  2.9O8OF-OI 
C3s  4.4704e  00 

DCk3s  1.0300E  01 
DK3s  1.Q300E  01 
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as  1.3000E  00 


Rlr  5.0000F-01 
Cl-  2.6000E  00 

OCKlr  2.2500E  00 
OKI:  2.2500F  00 


R2(l)s  7.1691E-01 
C2(l)=  1.8571E  00 

DCK2=  1.0525E  01 
DK2r  1.0525E  01 


r3x  2.9080E-01 
c3«  4.47o4p  00 

nC«3»  1.0300E  01 
DK3:  1.Q300E  01 


A=  1.300DE  00 


A:  1.3000E  00 


Rl=  1. 000OE  08 
Cl-  1.3000E-08 

DCKls  2.25QQE  00 
DKls  2.25Q0E  00 

Rls  1.0000E  08 
C^r  1.3000E-08 

DCKls  2.2500E  00 
DKl=  2.2500E  00 


R2(2)=  4.1542E  07 
C2(2)S  3.2o49E-08 

DCK2=  1.0525g  01 
DK2z  1.05?5e  01 

r2(1):  1.4338E  oe 
C2(l)=  9.2854E-09 

DCK2S  1.0525R  01 
DK2=  1.0525E  01 


»3E  5.8160E  07 
c3«  2.2352F-08 

nCK3=  1.0300E  01 
DK3s  1.0300E  01 

r3s  5.8160P  07 
c3«  2.2352p-08 

Dck3=  1.0300E  01 
nK3=  1.03Q0F  01 


a:  1.3500E  00 


4=  1.3500b  00 


Rl=  2.0000E  05 
Cl:  6.7500E-06 

DCKls  2.2500E  00 
DKl=  2.2500E  00 

Rl=  2.0000E  05 
Clz  6.7500E-06 

DCKlr  2.2500E  00 
DKls  2.2500E  00 


r2(2)S  8.2537E  0* 
C2(2);  1.495gE-05 

DCK2S  1.0525E  01 
DK2=  1.0525E  01 

R2(l):  2.5o52e  05 
C2(l)=  4.92g0E-06 

DCK2=  1.0525e  01 
DK2i  1.0525E  01 


r3s  1.1305F  05 
C3s  l.l941p-05 

nCK3=  1.0300E  01 
nK3=  1.0300F  01 

r3=  1.1305F  05 
C3=  l.l94ip_05 

DCK3r  1.0300F  01 
DK3s  1.Q300E  01 


A  =  1 .  35Q0E  .  00 


a:  1 . 3500E  on 


Rl=  2.5000E  05 
Cls  5.4000E-06 

DCKls  2.2500=  00 
PKls  2.2500E  00 

sir  2.5000F  05 
Clr  5.4Q00E-06 

DCKls  2.2500c  00 
DKls  2.2500E  00 


R2(2):  1.o317E  05 
C2(?)=  l.l966E-05 

DCK2:  1.0525E  01 
DK2=  1.0525E  01 

R2(l)s  3.1315b  05 
C2(l)=  3.9424E-06 

DCK2r  1.0525E  01 
DK2z  1.0525E  (11 


r35  1.4i32F  05 
c3:  9.553lp-06 

DCK3r  1.0300E  01 
DK3s  1.0300E  01 

r3-  1.4l32p  05 
C$  =    9.553lp-06 

ncK3r  1.0300E  01 
nK3=  1.0300E  01 


As  1.3500fc  00 


as  1.35Q0F  00 


fir  1.00'OOE  06 

Cls  1.3500E-06 

DCKls  2.2500E  Ot 
DKls  2.2500E  00 

Rl=  1.0000F  06 
Cls  1 . 350HF-06 

DCKls  2.250UE  00 
DKls  2.2500E  00 


r2(2)s  4.1268E  05 
C2(2)=  2.99l6E-06 

DCK2=  1.0525E  01 
DK2z  1.0525E  01 

p2(l)=  1.2526E  06 
C2(l)r  9.8560F-07 

DCK2=  1.0c525e  01 
DK2:  l.o525p  01 


93»  5.6526F  05 
c3r  2.3883F-n6 

DCK3=  1.0300E  01 
nK3r  1.0300E  01 

R3r  5.6526b  05 
c3:  2.3fl83F-06 

DCK3r  1.0300E  01 
DK3:  1.0 -''OOf  01 


t.-    1.35 OOP  00 


-lr  1  .5000t  05 
rlr  o.nnnnu-n6 


r2(2):  6.1903p  04 
r2(2)=  1 . oo44f-05 


r3=  8.4789F  04 
r3s  1.5o22c-05 
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A*  1.35Q0E  00 


DCKl*  2.2500g  00 
DKls  2.P500E  00 

Rls  1.5000F  0' 
Cl=  9.0000P-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 


DCK?*  1.0525g  01 
^OK2s  1.0525E  01 

R2(l)=  1.8789E  05 
C2(l)«  6.5707E-06 

DCK2a  1.0525E  01 
_ftK2a  1.0525E  01 


DCK3»  1.0300E  01 
{1*3.8.  1.030  0E  01 

r3i  8.4789E  04 
c3i  1.5922f-05 

DC*3i  1.0300E  01 
D*3S  1.0300E  01 


A«  1.3500E  00 


A:  1.3500E  00 


Rl=  S.0000E-01 
Qt*    2.7000E  00 

DCKls  2.2300E  00 
DKls  2.25Q0E  00 

His  5.0000E-01 
Cii  2.7000E  00 

DCKls  2.2500E  00 
DKls  2.2500E  00 


R2(2)i  2.0634E-01 
C2(2)i  5.9831E  00 


DCK2.  1. 05256  ol 
DK2b  1.0525E  01 

R2(l)t  6.2630E-01 
C2(l)»  1.9712E  00 

DCK2»  1.0525e  01 
DK2:  1.0525E  01 


r3i  2,s263f.0l 
C3«  4.7765E  00 

DC»<3i  1.0300E  01 
dk3«  1.Q300E  01 

p3,  2.8263e.01 
C3i  4.7765E  00 

OCK3i  1.0300E  01 
Dk3=  1.0300E  01 


As  1.3500P  00 


A:  1.35Q0E  00 


Rls  1.0000E  08 
Cl=  1.3500E-08 

DCKls  2.25QQE  00 
DKls  2.3500E  00 

Rls  1.0000E  08 
Clt  1.3500E-08 

DCKls  2.2500E  00 
DKls  2,?500E  00 


R2<2)»  4.1268E  07 
C2(2)s  2.9916E-08 

DCK2i  1.0525E  01 
DK2s  1.0525E  01 

r2<1).  1.2526E  08 
C2(l)t  9.8560E-09 

DCK2«  1.0525E  01 
DK2s  1.0525E  01 


r3«  5.6526F  07 
C3.  2.3ee3E-08 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 

r3«  5.6526E  07 
c3«  2.3ee3F-08 

nCK3i  1.0300E  01 
DK3s  1.03Q0E  01 


As  1.4000E  00 


As  1.4000E  00 


Rls  2.0000E  05 
Cls  7.0000E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 

Rls  2.0000E  05 
Cta  7.0000E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 


R2(P)s  8.?722E  04 
C2(2)»  1.3877E-05 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 

R2(l)i  2.1940E  05 
C2(l)«  5.2324E-06 

DCK2i  1.0525E  01 
DK2t  1.Q525E  01 


r3«  1.1067E  05 
C3»  1.265tE-05 

DCK3i  1.0300E  01 
DK3s  1.03Q0E  01 

r3»  1.1067E  05 
C3i  1.2651F-05 

DCK3«  1.0300F  01 
DK3s  1.Q300E  01 


As  1.4000E  00 


as  1.4000E  00 


Rls  2.5000E  05 
Cl=  5.6000E-06 

DCKls  2.25Q0E  00 
DKls  2.2500E  00 

wis  2.5000E  05 
Cls  5.6000E-06 

DCKls  2.2500E  00 
DKls  2.2500E  00 


R2(2)s  1.0340E  05 
C2(2)s  1.1102E-05 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 

R2(l)s  2.7424E  o5 
C2<1)»  4.i859E-06 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 


r3i  1.3g33E  05 
C3i  1.0120E-05 

DCK3s  1.0300E  01 
DK3s  1 .  0300F  01 

r3.  l,3833p  q5 
c3t  1.0120E-05 

DCK3s  1.0300E  01 
DK3s  1.0300F  01 
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A=  1.4000E  00 


A:  1.4000E  00 


Rli  1.0000E  06 
Cla  1.4000F-0* 

DCKls  2.2500E  00 
DKl=  2.2500E  00 

Rls  1.0000E  06 
Cls  l.«O0OE"0* 

DCKl=  2.2500E  00 
HKls  2.2500E  00 


R2(2)r  4-.1361E  05 
C2(2)s  2.7755E-06 

DCK2=  1.0525g  01 
D«2s  1.0525E  01 

R2(l)s    1.0970E    06 
C2(l)r    1.04>65E-06 

DCK2:  1.0525e  01 
DK2=  1.0525E  01 


p3.  5.5334E  o5 
C3S  2.53olF-06 

DCk3j  1.0300E  01 
DK3s  1.0300E  01 

R3«  5.5334E  05 
p3i  2.5301F-06 

DCK3=  1.0300F  01 
nK3=  1.0300E  01 


A=  1.4Q0OE  00 


A=  1. 4000F  00 


Bl:  1.5000E  05 
Cls  9.3333g-o6 

DCKls  2.2500=  00 
DKl=  2.2500E  00 

Rls  1  .5000E  05 
Cls  9.3333g-o6 

DCKls  2.2500F  00 
DKl=  2.?500E  00 


r2(2)s  6.2q4ic  q4 
C2(2)r  1.85o3e-o5 

DCK2=  l.o5?5g  01 
DK2=  1.0525E  01 

R2(l)s  1.6455E  05 
C2(l)a  6.9765E-06 

DCK2:  1.0525F  01 
Dk2=  1.0525E  01 


r3«  8.300lF  04 
C3=  1.6867F-05 

DCk3=  1.0300E  01 
OK':  1.Q300E  01 

r3s  8. 3 00  IF  04 
C3s  1.6867F-05 

DCK3b  1.0300E  01 
[)K3t  1.0300F  01 


As  1.4000E  00 


»=  1.  4pp0E  00 


pis  5. 000CP-01 
Cl=  2.80006  00 

DCKls  2.2500=  00 
DKls  2.25P0F  00 

wis  5  .  0000F-01 
C-ls  2.3000F  00 

DCKls  2.25(10  =  00 
HKl:  ?.?5nOF  00 


R2(2)=  ?.0680F-01 
C2(2)r  5.5509E  00 

DCK2;  1.0525F  01 
DK2=  1.05P5E  01 

r2(1)=  5.4fl49t-01 
C2(l>=  2.D930E  00 

DCK2=  1.0525E  01 
DK2r  1.0525E  01 


r3=  2.7667F-01 
c3s  5.0602p  00 

DCK3=  1.0300E  01 
DK3r  1.0300E  01 

r3=  2.7667F-01 
c3r  5.0602F  00 

DCK3=  1.0300E  01 
D«3=  1.Q3Q0E  01 


A=  1.4  0  OOF  00 


A=  1.4000F  00 


wis  1.0000E  08 
Cls  1.4Q00F-08 

DCKls  2.250HE  00 
OKls  ?.?500F  00 

Kir  1. 0000F  08 
Cls  1.4Q00F-08 

DCKls  2.2500=  00 
DKls  2.P500F  00 


r2(2)=  4.1361E  07 
C2(2)s  2.7755F-08 

DCK2:  1.0525E  01 
DK2=  1.0525E  01 

r2(1)S  1.0970E  08 
C2(l)=  1.0465E-08 

DCK2i  1.0525F  01 
DK2i  1.0525E  01 


p3s  5.5334c  07 
c3=  2.53plE-08 

DCK3r  1.0300E  01 
DK3=  1.0300F  01 

r3s  5.5334c  07 
c3=  2.5301F-08 

Dck3s  1.0300E  01 
DK3-  1.0300E  01 


as  1.4500F  00 


i=  1.4500E  On 


wl=  2  .  0000E  05 
Cls  7.2500E-06 

DCKls  2.2500E  00 
pKlr  ^.2500E  00 

ols  2. 000  OF  05 

H:  7  25nnc-n6 


R2(2)s  a.  3s55E  04 
C2(2)s  1.2762E-05 

DCK2s  1.0525E  01 
D«2=  1.0525E  01 

R2f 1 ) =  1.Q180E  05 

r2(  l )  z    5  .  57QISC-  n* 


r3=  1.0896F  05 
c3=  1.3308F-05 

Dck3=  1.0300F  01 
n«3r  1.Q300F  01 

°3r  1 . 0896E  05 
r  3  s  1.3  3n«c-05 
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DCKls    2.2<5nOE    00 
DKl:    2.25006    00 


DCK2t  1.0525E  01 
_DK2  =  1.05?5E  01 


DCK3s  1  .  0300F  01 
DK3:  1.Q300E  01 


A=  1. 45Q0E  08 


A=  1.45Q0F  00 


ol:  2.5000E  05 
Cl=  5.8000F-06 

DCKl  =  2.2500E  00 
DKl:  2.25006  00 

Bis  2.5000E  05 
Cl=  5.e000F-06 

DCKls  2.2500E  00 
DKl=  2.2500F  00 


r2(2,=  1.0482E  05 
C2(2):  1.0210E-05 

DCK2r  1.0525E  01 
DK2=  1.0525E  01 

r2(1)-  2.3975E  05 
C2(l)=  4.4636E-06 

DCK2i  1.0525E  01 
DK2r  1.0525E  01 


r3e  1.362op  05 
c3s  1.0646F-05 

DCK3s  1.0300E  01 
nK3=  1.0300E  01 

r3s  1.3620F  05 
C3s  1.06"6F-05 

DCK3:  1.0300F  01 
DK3r  1.Q300E  01 


a=  1.45Q0E  00 


a=  1.45Q0E  00 


His  1. 0000F  06 
Cls  1.45006-0* 

DCKls  2.2500E  00 
DKls  2.25006  00 

pl=  1.0000F  06 
Cls  1. 45006-06 

DCKls  2.25Q0E  00 
DKl:  2.2500E  00 


R2(2)=  4.1927E  05 
C2(2)=  2.5524E-06 

DCK2=  1.0525E  01 
DK2=  1.05256  01 

r2(1)=  9.59006  05 
C2(l)=  1.11596-06 

DCK2s  1.0525e  01 
DK2=  1.q525e  01 


r3«  5.4480F  05 
C3s  2.6615F-06 

nCK3=  1.0300E  01 
nK3=  1.0300E  01 

r3b  5.4480E  05 
c3s  2.6615f-06 

DCK3s  1.03006  01 
0K3i  1.Q300E  01 


As  1.45006  00 


A=  1.4500E  00 


Rls  1.5000E  05 
Cls  9.6667E-06 

DCKls  2.2500E  00 
IlKls  2.25006  00 

pis  1 .5000E  05 
Cls  9.6667E-06 

DCKls  2.2500E  00 
DKl:  2.25006  00 


R2(2)=  6.289lE  04 
C2(2)=  1.70166-05 

DCK2=  1.05256  01 
DK2s  1.0525E  01 

R2(l)r  1.43856  05 
C2(l)=  7.4394E-06 

DCK2:  1.05256  01 
DK2:  1.0525E  01 


p3i  8.1720F  04 
C3s  1.7744E-05 

DCK3s  1.0300E  01 
DK3s  1.03006  01 

r3s  8.17206  04 
C3s  1.7744E-05 

DCK3s  1.03006  01 
DK3:  1.03006  01 


as  1.45006  00 


as  1.4500F  00 


ol:  5.0000F-01 
Cls  2.9000F  00 

DCKls  2.25006  00 
DKl:  2.  ?5fl0t  00 

ols  5. OOOHF-Ol 
Cl  :  2.90006  00 

DCKls  2.2500=  00 
DKl:  2.P5006  00 


r2(2)=  2.0964e-01 
C2(2)r  5.1Q486  00 

DCK2:  1.05?5c  01 
DK2:  1.05256  01 

r2(1)=  4.795oe-01 
C2(l)s  2.23186  00 

DCK2:  1.0525e  01 
DK2:  1.05256  01 


r3s  2.7240F-01 
f-3:  5.3231F  00 

DCK3:  1.03006  01 
nK3:  1.03006  01 

«3r  2.7240F-01 

c3:  5.3231F  00 

nCK-?:  1.0300F  03 
n<3:  1.030DF  01 


As  1.45006  00 


Rls  1.00006  08 
Cls  1.45006-08 

nr<ls  2.25noc  no 


B2(2)s  4.  19276  07 
C2(2):  2.55?4E-08 

nrK2;  1.05?5p  nl 


ois  5.44806  07 
c3s  2.66156-08 

nrK3s  1.03P0F  01 
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OKlr  2 . ?50  OP  00 


DK2:  1.0525E  01 


0*3=  1.Q300E  01 


a=  i.«5qoe  on 


si:  l.OOOOE  08 
Clr  1.4500E-08 

ncKl=  2.2500E  00 
nKls  2.2500E  00 


R2(l)=  9.5<}O0E  07 
C2(l)r  1.1139E-0S 

DCK2s  1.0525E  ol 
DK2=  1.05?5e  01 


r3»  5.4480e  07 
c33  2.6615E-08 

DCK3r  1.0300E  01 
pK3r  1.Q300E  01 


A=  1.5000E  00 


As  1.5000E  00 


Rl:  2.0000P  05 
Cl=  7.5000E-06 

DCKls  2.2500E  00 
DKl=  2.2500E  00 

Rls  2. 0O0OF  05 
Cl=  7.5000E-06 

DCKls  2.2500E  00 
PKls  2.?5ooE  00 


R2(2)r  8.6443E  0* 
C2(2)s  1.1568E-05 

DCK2r  1.0525g  01 
DK2=  1.Q525E  01 

R2(l)=  1.6625E  05 
C2(l)s  6.0150E-06 

DCK2=  1.0523e  01 
DK2=  1.Q525E  01 


r3*  1.0778F  05 
C3.  1.39l7F-05 

nCK3s  1.0300E  01 
nK3=  1.0300E  01 

r3=  1.0778F  05 
C3»  1.3917F-05 

DCk3s  1.0300E  01 
nK3s  1.o3qoe  01 


A=  1 .5000F  00 


A=  1. 5000E  00 


Jl:  2 .5000E  05 
Cl=  6. OOOOE-06 

DCKls  2.2500E  00 
DKls  2.?5noF  00 

Bis  2.5000F  05 
Cls  6.0000E-06 

DCKls  2.2500E  00 
nKls  2.25onE  00 


R2(2)=  1.C905F  05 
C2(2)=  9.2547F-0* 

DCK2r  1.0525F  ol 
DK2=  1.o5?5e  01 

R2(l)=  2. 0781E  05 
C2(l)s  4.fli20F-06 

DCK2r  1.0525E  ol 
DK2r  1.0525E  01 


r3=  1.3473e  o5 
C3=  1.1133f-05 

DC«3s  1.0300F  01 
nK3=  1.0300E  01 

r3s  1.3473e  o5 
c3s  l.ll33F-o5 

DCK3s  1.0300F  01 
nK3r  1.03Q0F  01 


A=  1.5000E  00 


A=  1 . 5  0  G  0  F  0  0 


Rl:  l.OOOOE  06 
Cls  1.5000F-0<S 

OCKls  2.2500E  00 
nKlr  2.2500E  00 

si:  l .  onnoF  06 

rl  -  1  .  5000E-06 

nCKls  2.2500E  00 
nKls  2.2500E  00 


r2(2):  4.3?2lE  05 
C2(2)=  2.1il37E-06 

DCK2s  1.0525E  01 
DK2z  1.0525E  01 

R2(l)=  8.3126E  05 
C2(l)=  1.2n30F-06 

DCK2s  1.0525E  01 
DK2=  1.0525E  01 


r3=  5.3b92f  05 
C3S  2.7833F-06 

dCk3=  1.0300E  01 
OK3=  1.0300E  01 

r3-  5.3892E  05 
c3:  2.7833F-06 

nCK3s  1.0300E  01 
nK3r  1.03Q0E  01 


as  1.5000R  00 


i:  1.5000F  00 


Rl:  1. 5000£  05 
Cls  1.0000F-015 

DCKls  2.2500E  00 
DKls  2.2500E  00 

si:  1  .  5000E  05 
Cls  1 . 0OnOE-05 

DCKls  2.25Q0E  00 

QKlr  2t?500E  on 


r2(2)=  6.4832e  0< 
C2(2)=  l.(5424E-05 

DCK2:  1.0525e  01 
DK2:  1.0525E  01 

r2{1)=  1.2469E  05 

C2( 1):  a.  0200F-06 

nCK2:  1.0525e  ol 
DK2=  1.05?5e  01 


r3s  8.0838F  0A 
c3:  1.8556F-05 

DCK?:  1.0300F  01 
nK3=  1.0300E  01 

R3=  8.0838F  04 
c3=  1.B556F-05 

DCK3s  1.0300F  01 
nK3=  1.0300c  01 


83 


A:  1.5000E  00 


A:  1.5000E  00 


rIs  5.0n0OE-0l 
Cli  3.0000P  on 

DCKl:  2.2500E  00 
DKls  2.2500E  00 

Rls  5  .  Cin00E-0l 
Cl=  3. OO0OE  00 

DCKl:  2.?'5on=  00 
DKls  ?.?5nni-  00 


r2<?)=  2.16HE-01 
C2(2)=  4.6273E  Of) 

DCK?=  1.0525E  01 
DK2=  1.0525E  01 

r2(1)=  4.1563E-01 
C2(l)r  2.4060E  00 

DCK2=  1.0525P  01 
DK2r  1.£)525e  01 


r3s  2.6946E-01 
c3i  5.5667E  00 

nCK3=  1.0300E  01 
DK3r  1.0300F  01 

r3s  2.6946E-01 
C3a  5.5667F  00 

nC*3  =  1.0300E  01 
nK3;  1.q3ooc  01 


A=  1.5000E  00 


As  1.5000E  00 


Bl=  1. 0000E  08 
Cl=  1.5000E-08 

DCKl:  2.2500E  Q0 
DKl=  2.2500E  00 

Rls  1.0000E  08 
Cls  1.5000E-08 

DCKl:  2.2500=  00 
DKls  2.2500E  00 


R2(2)=  4.3221E  07 
C2(2)s  2.3137E-08 

DCK2s  1.0525E  01 
DK2s  1.0525E  01 

r2(1)=  8.3126E  07 
C2(l)s  1.2o3oe-08 

DCK2s  1.0525£  01 
DK2:  1.0525F  01 


r3i  5.3892e  07 
c33  2.7833f-08 

DCK3s  1.0300E  01 
nK3:  1.0300F  01 

r3s  5.3a92p  07 
c3=  2.7833c-08 

HCK3s  1.0300E  01 
DK3s  1.0300F  01 


As  1.5500F  00 


4=  1.55Q0E  00 


Pi:  2. 0000E  05 
Cls  7.75Q0E-06 

DCKls  2.?500E  00 
DKls  2.2500E  00 

Rls  2.0000E  05 
Cls  7  .  7500E-06 

DCKl:  2.2500E  00 
DKls  ?.?500E  00 


r2(2)=  9.2118F  04 
C2(2)s  1.0167E-05 

DCK2s  1.0525e  01 
DK2s  1.0525E  01 

R2(l)s  1.404lE  05 
C?(l)s  6.6698E-06 

DCK2s  1.0'525f  01 
DK2=  l.n525E  01 


r3«  1.070<F  05 
C3s  1.4481E-05 

DCK3s  1.0300E  01 
DK3s  1.0300E  01 

r3s  1.0704F  05 
C3»  1.44B1F-05 

nCK3s  1.0300E  01 
DK3s  1.q3qof  01 


A:  1.55Q0E  00 


a:  1.5500F  00 


Rls  2. 5000F  05 
Cls  6.2000F-06 

DCKl:  2.2500E  00 
DKls  2.P500F  00 

Rls  2.5nnOE  05 
Cls  6.2000F-06 

DCKl:  2.25005  00 
DKls  2.2500E  00 


R2(2)s  1.1515E  05 
C2(2)s  8.1332F-06 

DCK2:  1.0525E  01 
DK2=  1.0525F  01 

R2(1)S  1.7551F  05 
C2(l)s  5.3359E-06 

DCK2s  1.0525F  01 
DK2s  1.05?5e  01 


r3s  1.3380E  05 
C3=  1.1585E-05 

DCk3=  1.0300F  01 
OK3s  1.0300F  01 

qi=    1.3380F  05 
C3r  1.1585F-05 

DCK3:  1.0300E  01 
nK3:  1.0300c  01 


As  1.5500E  00 


as  1.55Q0E  00 


pis  1.0000=  06 
Cls  1.5500E-06 

DCKl:  2.2500E  00 
DKls  2.?5ooc  00 

els  1. 0000F  0* 
Cl=  1.55Q0E-06 

ni-Kis    2  .  ?5nn=    nn 


R2(2)s    4.6059e    05 
C2(2)s    2.0333F-06 

DCK2s    1.05?5e    01 
DK2s    1.0525E    01 

R2(l)=    7.0206E    05 
C2(l)s     1.3  34qf-06 

ncK2:  1  .o525p  nl 


r3s  5.3518F  05 
C3s  2.8962F-06 

DCK3s  1.0300F  01 
nK3:  1.0300E  01 

r3:  5.3518E  05 
C3s  2.8962F-06 

nrK3=  1  .  o3on=  01 
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DKlr  2.2500E  00 


DK2:  1.05?5e  01 


DK'=  1 .  0300F  01 


A»  1.5500E  00 


A  =  1.55Q0F  00 


Rls  1.5000F  05 
Clr  1.0333e-05 

DCKl:  2.2500E  00 
DKlr  2.2500E  00 

Rls  1 . 5 n  0  0 E  0  5 
Cls  1.0333E-05 

DCKl=  2.25Q0E  00 
DKls  2.?500E  00 


R2(2)«  6.9089E  04 
C2(2)r  1,355'5e-05 

DCK2i  1.0525E  01 
DK2=  1.05?5e  01 

R2(l)r  l.o53lg  05 
C2<1>*  8.893lE-06 

DCK2=  1.0'525e  01 
DK2:  1.o5?5e  01 


r3«  8.0278E  04 
c3«  1.9308E-05 

DCK3=  1.0300E  01 
qk3=  1.0300E  01 

r3«  8.0278F  0* 
c3«  1.9308E^05 

nCK3=  1.0300E  01 
nK3r  1.0300E  01 


A=  1.55Q0F  00 


A:  1.5500E  00 


Rlr  5. OOOOfc-Ol 
Cl=  3. 1000E  00 

nCKl=  2.25r)0E  00 
DKl=  2.?500E  00 

olr  5. 0000E-01 
Clr.  3.1000E  00 

DCKlr  2.2500E  00 
DKlr  2.2500E  00 


r2(?)=  ?.3o3oe-01 
C2(2)=  4.0666E  00 

DCK2r  1.0525E  01 
DK2r  1.0525E  01 

R2(l)=  3.5103f-01 
C2(l)r  2.6679E  00 

DCK2r  1.05?5f  01 
DK2z  1.05?5f  01 


p3r  2.6759F-01 
c3=  5.7924p  00 

DCK3r  1.0300E  01 
nK3r  1.0300E  01 

«3r  2.6759F-01 
C3=  5.7924F  00 

DCK3r  1.0300F  01 
nK3r  1.Q300E  01 


Ar  J..55Q0E  00 


Ar  1.55Q0F  00 


R 1 r  1 . 0  0  0  0  E  08 
Clr  1.5500E-08 

DCKl=  2.2500=  00 
DKlr  2.2500E  00 

uls  1 . 0000F  08 
Clr  1.5500F-08 

nc<l=  2.2500E  00 
DKlr  2.2500E  00 


R2(2)r  4.6()59E  07 
C2(?)r  2.0333E-08 

DCK2r  1.0525c  oi 
DK2r  1.Q525E  01 

R2(l)s  7.  0206F  07 
C2(l)r  1.3340E-0B 

DCK2r  1.0525E  01 
DK2=  1.0525E  01 


r3=  5.3518F  07 
c3r  2.8962F-08 

DCK3=  1.0300E  01 
DK3r  1.0300F  01 

r3=  5.3518E  0  7 

c3=  2.8962fc-08 

DC<3r  1.0300E  01 
DK3:  1.Q300E  01 


8<r 
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PROGNA*   DRAWCCMP 

DIMFNSION  DRAW^1(900).DRAwp21(900),DRAWR22<900>, 
lDRAWP'(9COi.D^AWCl(9f'0).DRAWc21(900).DRAWC22(900),DRAWC3(900)1 

2nwAWA<900) , n i Tue ( 12  > 

DO    999    J=l,900 
DRAWRM  JJsO.O 

DNAWR21 ( J):0 . 0 
DRAWR2?< J)r0  .  0 
DRAWR3( j)s0  .  0 
DRAHCl ( J)  rO  .  0 
DRAWC21( J)  =0.0 
DRAWC22( J)rO.O 
DWAWC3( j)  :0  .  0 

099  DRAWA( J)  =  0  .  0 

read  5o,dki.d<?.dk3.starta.c;tepa.stopa 

50  F0WMAT(6Fln ) 

READ  60, Rl 
60  FOHMAT(FlO) 

A  :  S  T  A  «~  1  A 

1  =  1 

100  YsDK?/ A-A- ( 2 . «CKl )/( A**2) 
IF( v )  101,101,102 

1  01  A  =  A*S^  t-PA 

IFCSTOPA-All?'',  125,100 

102  XsDK3-?.«a-0K1/'A**2) 
I F  c  X )  101,101,103 

103  Tl3=x-Y 

IFl Tl3)101, 101,  104 
10"  W= I V**2)-(4 . *QKl*rl3) /( A**2) 
IF(W>  101,105,  105 

105  CUNT  I NUE 
C3-T13/P1 
Cl=A/Pl 

R21=  (  Y  +  SORIf"  (  A)   W(2,*C31 
C2l  =  DKl/( ( A**2)*R2l  i 
R22:  (  y-SORti.-  (  ^  )  )  /  (  ?  .  „o3  ) 

106  c22=DKl/( ( a*.2>»^22) 
nKAk.Pl  (  I  ) sPl 
PRAWR21 ( I ):R2l 
DkAKR?'(  I  )  z-i?-i 

D  B  A  W  P  ?  (  I ) =  R  3 
D*AkC1  (  I >  =C 1 

riBAwc2i  ( i )  -c?i 

DRAKC2?( I )rc22 

nwAujc^f  i  )  =  c3 

PPAWA ( I )      :t 

I --1*1 

GO  TO  101 
125  CONTINUE 

N  U  f  -  6  9  P 

DO  7^0  JK=1.12 
700  IT 1TLE' J*  )  =  8* 

I  TITLE!  1)=  8"'_  ILL  IS. 

!  TITLff  <2)=  a^  .  .  A. 

I T 1 T  L  h ( 3  )  s  3h  thESI^ 
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ITITLE<4)r    9H         PLOT 
ITITLF(5)r     8H0F     P2(l) 
IT  I TLE«  6) *    8H.     R?(2) . 

ITITLF(7)r     8H     AND    P3 

1TITLE<8>=«HVS     UME 

!TITLE(9)s  8HCCNSTANT 

I  TITLE( 10 ) =8H.  A.  FOW 

ITITLE(H)=8H  CON5TAN 

ITITLE(12)=8HT  Rl 

LABEL=4HR2l 

CALL  D«  AW(NUM.CRA*/A.DRAW^2l. 1,0, LABEL. 1UTLE, 0,0, 0,0, 
in, 0,0, 0,0, LAST) 

LABEL=4HRl 

CALL  DRAW(NUM.r><A^A.DRAW>Jl, 2,0. LABEL. ITITLE, 0,0, 0,0, 
10,0,0,0,0,1  AST) 

LABEL=4HR22 

CALL  DWAW(NIJm.cRa^A.D:?AWR22,  2,0, LABEL- I  TITLE.  0,0,  0,0, 
10.0,0.0, 0,1  AST) 

LABEL=4HR3 

CALL  DWAW(NiJM,r^A^A.DRAW»3,3,n.LA8EL,  J  T  I  TLE ,0,0,0,0, 
10,0,0.0, O.LAS  r ) 

PO  701  JKs5,l2 
701  ITITLE( JK)r8H 

I  T  I  TLE  <  5 ) =  SHOP  Cl.C2 

ITITLE(6)s  8H(1 ) , 

!TlTLE(7):  8WC2(2) ,  a 

I  T  I  TLE ( 8 ) =  8H\r  C3  VS 

IT  I  TLF(9)=  8H  T If C  CC 

IT  l  TLF(  10  )  =8W\<:TaNT  . 

ITiTLE(H)r8HA,  F^R  C 

I  T  ITLE(l2)::<}l-OST  .  Wl 

LABEL=4HC2? 

CALL  DPAW(N'JM,r-»A^A.D3AniC2?, 1,0, LABEL.  I  TITLE.  0  .  0.0,0, 
10,0,0,0,0,1  ASf> 

LABbL=*HCl 

CALL  D^AW(MJM,C»A^A.DRAWCl, 2.0, LABEL, I T I TLE ,0,0,0.0, 
ir.n,0,0,0,LA^T) 

LABEL=4HC21 

CALL  DMAW(NiJM,rRAWA.DRA^c2l. 2,0, LABEL,  I  TITLE. 0,0, 0.0. 
10,0,0,  0,0, LAST! 

LaBEL=*HC3 

CALL  DKAWtK'J^.CPA^A.D^A^C3. 3.0. .ABEL,  I  T  I  TLE  ,0.0.0.0. 
in,  n,  0,0,  0,1  AST) 

FND 
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APPENDIX  IV 


CONTINUANT  EXPANSIONS 


It  was  previously  stated  that  the  expansion  of  a  continuant  can 
be  formulated  as  the  evaluation  of  a  signal-flow  graph  representation 
based  on  Euler's  rule.   The  method  is  especially  useful  when  working 
with  RC  ladder  networks  since  all  the  b  coefficients  of  the  continu- 
ants are  unity.   The  signal-flow  graph  for  the  network  would  be  as 
shown  in  Fig.  11. 


Fig.  11.   Flow  Graph  of  a  Continuant 
Rules  for  evaluation  of  the  signal-flow  graph  are  as  follows: 

1)  The  first  term  is  given  by  the  product  of  all  adjacent  R  and 
sC  transmissions. 

2)  The  second  term  is  given  by  the  sum  of  products  formed  from 
the  first  term  by  omitting  adjacent  pairs  of  R  and  sC  trans- 
missions, taken  one  pair  at  a  time,  in  every  possible  way. 

(This  is  equivalent  to  setting  the  transmission  of  the 
omitted  pair  equal  to  unity.) 

3)  The  third  term  is  given  by  the  sum  of  products  formed  from 
the  first  term  by  omitting  adjacent  pairs  of  R  and  sC  trans- 
missions, taken  two  pairs  at  a  time,  in  every  possible  way. 
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4)   The  nth  term  is  given  by  the  sum  of  products  formed  from 
the  first  term  by  omitting  adjacent  pairs  of  R  and  sC 
transmissions,  taken  (n-1)  pairs  at  a  time,  in  every  pos= 
sible  way. 
In  forming  the  sums  of  products  specified  in  the  above  four  rules, 
it  is  convenient  to  think  of  the  omissions  of  a  transmission  of  an 
adjacent  R  and  sC  pair  as  being  equivalent  to  "bridging  around1"  this 
pair  with  a  path  of  unity  transmission. 
The  continuant  expansions  of 
<y.)   -  k(R1,sCn) 

for  values  of  m  and  n  up  to  m  =  8  and  n  =  4  are  as  follows? 

Q0(s)  =  1 

Q^s)  =  Rt 

Q2(s)  =  (R1C1)s  +  1 

Q3(s)  =  (R1C1R2)s  +  (Rx+  R2) 

Q,(s)   =   (R1C1R2G2)s2  +  (R1C1+  R^^  R2C2)s  +  1 

2 

Qc(s)      =      (R1C1R0C,RQ)s     +      (R1C1R0+R1C1R,+R1CoRQ+R0C0R,)s 
5  112/3  112      113      123      223 

+    CR^R^     R3) 

Q6(s)   =   (R1C1R2C2R3C3)s3  4   (R^R^^  R1C1R2C3+  R^R^ 

2 
+  R  C2R3C3+  R2C2R3C3)s   (R.CL+  RjC2+  R^C^  R2C'2 

+  R2C3+  R3C3)s  +  1 
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Q?(s)   =   (R1C1R2C2R3C3R4)s   +   (R1C1R2C2R3+  R  C^C^ 


+  R1C1R2C3R4+  R1C1R3C3R4+  R-^R^R^  R2C2R3C3R4)s 


+  (R1C1R2+  R1C][R3+  RjCjR^I-  R^^+R  C2R4+  RiC-jR, 


+R2C2R3+  R2C2R4+  R2C3R3+  R-C-R^s   +   (R  +  R  +  R  +  R.) 


Qg(s)   =   (R1C1R2C2R3C3R4C4)s   +   (R1C1R2C2R3C3+  R^R^R^ 


"   Rn  C,  R0C~R/  C. t   R,  C.  Rr,C« R/ C.  "t   R- Cn  R0C0R,  C, 
112244        112344        113344 


+  R1C2R3C3R4C4+  R2C2R3C3R4C4)s        +      (R1C1R2C2 


t   Rn  C.  R0C-t   R.C»R0C/T   R,  C.RoC^t   RtC-R^C/T   R.  Cn  R,  C. 
1123        1124        1133        1134        1144 


■   Rn  C0    R/>C/,t   R.C0R0C/T   RnC0R/C/  i    R-.C-.R.  C#t   RoCoR0C0 
12,3.  3        1234        1244        1344        2233 


+  R2C2R3C4+  R2C2R4C4+  R2C3R4C4+  R3C3R4C4)s 


+   (Rj^C-H-  R,C2+  RXC3+  RXC4+  R2C2+  R2C3+  R2C4+  R3C3 


+     R0C,+  R,C,)s  +  1 
3   4        4   4 
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